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AhsI  met 

Ihc  rnliusi  II.  ss  of  iMo\  iii^-l>aiik  nnilli|>l(;-  iiioili-l  a.la|iliv.  i  sl  iiuat  i.  .n  (  M  M  A  I  ' )  .aii'i  .  .  :iii  .  .1 

iMMAC)  algorithms  is  an.alyzc.l  in  this  lln-sis.  Tin'  misinalcli  of  a  Al-siali  iiiiili  m...|.[  mi  i  a 

I'.-siati  liltiT  mo.li'l  ari'  evaltialt'd  on  tin'  liasis  of  M  M  A  1 1/M  M  A<  '  porforman.-" 

A  nunli'l  il.'vrlop.al  using  finil"  I'lenicnt  analysis  is  usad  lo  approxiinali'  a  lai  spa.-r  si  i an  i  m  . 
will. 'll  has  a  hugi-  rontral  huh  with  aplx'tnlagns  radiating  out  from  it.  1  In'  mass  ,  •!  ih.’  Iml.  is 

.■oiisidurod  to  ho  much  largi'r  than  tlio  mass  of  tin.'  Ih'xil'l''  st  nni  nr.  .  lli.*  m.'.l.l  is  1 

III  physical  coor.liiiali's  and  then  l  ransr()r"n-d  intci  in...lal  coordina!. .  !..  .Tiam  a  i'..iii..  .1  i.i'  r 

lilt,  r  model ,  tin'  liwt  In  id  of  singular  perl  nrliat  ions  is  used .  1  In-  a.-i  nal  p.  .'i  i  n  .iis  and  \  i  n  -  .  .1 

'.  :ii'i.  .IIS  physi.'al  points  on  I  he  st  rnct  nre  are  n.sed  in  tin  -  e  vain  at  n.n  ■  .1  i  In  M  .\1  .\  li  /  M  \1  A(  '  ihg-i  ii  hm 
perlVirmance. 

1  ho  rosnlls  ol  this  study  of  model  mismatching  indicates  that  the  MMAli  pr.wi.les  ac.  iirato 
position  and  velocity  t'si.imtttcs  ('ven  in  the  f;ice  of  a  ti-stal.'  to  'Jl-.siate  m..,!.  1  mi.smat.  li.  When  a 
iioii-adcpt  ivo  nit.er  is  tised  with  'i  mismatcln'd  parameter  location,  the  performance  sulleis  slightU'. 
I  he  use  of  an  adaptive  estimator  iloes  |>r(iviih'  imprcived  performance  in  the  fa.'e  of  uin'.  ilain 
|.aramefer  location.  Stahle  control  was  ohiained  with  the  use  of  M.M.\('.  lor  the  case  of  noti- 
a.laplive  filt.'r  ami  mismtitched  ptirtmieli'r  location,  the  control  algorithm  l.ehaved  in  a  possihly 
destructive  manner.  By  allowing  the  lilli'r  lo  tidapi  to  the  initial  |.arame|i'r  location,  the  ,M.M,\(' 
algorit  Imi  [iroviih-.l  st  aide  cont  rol  of  t  he  si  rucl.iiri'.  even  m  t  he  fare  of  large  disi  nrhances. 


HOlH’S  l  NKSS  OF  A  MOVIXO-BANK 
.Ml  i;riBLF  MODFL  ADAPTIVF  COX^FHOLLFH 
FOB  A  I.ABC.’F  SPACF  STBfCTFKF 

I.  Inirodiiclion 

A  ,sil ii;il ion  rxists  in  nst iiiint ion  ;iiul  control  design  prohli'ins  where  tin'  |niraii h  I ms  descrilnng 
llie  system  are  not  ki\own  witli  certainty.  Additionally,  the  i)arani('t('r.s  may  Ijc  varying  as  a  rnnetion 
of  t  ime,  such  as  duo  to  t  lie  (h'plotion  of  a  fuel  cell  or  t  he  jet  I  i.soning  of  external  st  ores  on  an  aircralt . 
•Anotlier  possibility  is  th;it  the  system  ptiranieti.'rs  may  undi'rgo  a  jump  change  (a  radical  change 
from  otie  value  to  a  not  her)  dtti'  to  a  host  ile  environment  alfect  ing  t  In'  ph\'sical  si  met  ure  or  st  luci  nral 
failure  due  to  an  itii|ierfertion  or  structural  fatigui'.  Due  to  these  iiossihh'  parameter  \alue  changes, 
adaptive  estiitiation  of  the  systetn  slates  and  adaptive  control  might  he  tise(l  in  controlling  tin' 
system's  perfonmmee.  Aitolher  [irolih'tn  is  that  using  a  complete  mathematical  (or  "tritth")  model 
would  create  an  iinheartihle  computational  load  dm?  to  the  liigli  order  of  tin'  model.  Therelori'.  one 
Usually  chooses  to  usi'  a  rediiced-ordi'r  model  for  tin?  basis  of  synthesizing  a  lilter  and/or  c(.int  roller, 
and  so  robustness  of  the  ada[)tatioti  process  to  unmodelled  effects  becomes  an  important  issue. 

This  tliesis  tises  Kaltnati  ftlters  for  the  estimation  of  the  system  states  of  a  hirge  llexilile  space 
structure.  The  technique  of  estimatitig  the  system  states  is  referred  to  as  Mtiltiple  Model  Adaptive 
Ksfimation  (MiMAE),  having  the  structtire  of  a  bank  of  Kalman  filters  (each  based  upon  a  single 
discrete  value  of  the  parameters)  and  then  forming  the  final  state  estitnate  as  a  probabilistically 
weighti'd  sum  of  the  individual  filter  state  estimates.  Tin'  conc('|it  of  a  "'movingdrank  '  .MMAf. 
algorithm  is  that  not  all  possible  elemental  Kalman  lilU'rs  are  maintained  in  the  bank  at  all  times, 
but  only  a  dyntmiically  redeclared  subset,  such  as  those  corresponditig  to  discri'te  parameter  values 
clnsrst  to  the  ctirrent  parameter  estimate.  The  estimated  states  are,  m  turn,  used  a.s  iii|)uts  to 


a  r(Mil  ri  )llfr.  This  roiUrol  jiroc'css  is  rcfi'rnHl  to  as  Miillipli'  M(i<li-1  Ailaplivr  ('(inlr;)!  (MM.K). 
I’ll.'  MMAF,/.MM.\( '  alg(uit Inns  (F'volopfxl  in  this  rosnarcli  aro  nscal  to  comIioI  a  large  Ilexilil.'  sjiaia- 
St  nnt  lire'.  I' he  algor  it  Inns  have  evolvi'd  over  a  live  year  period  with  research  perforined  hy  llnii  /.  [.Sj . 
Filios  [it],  Karnick  [11],  Lashh'e  [15],  timl  \'an  Her  Werkeii  [2l|.  In  [tarticular,  the  work  |m  iliiriii' d 
ill  tin'  last  thi'sis  I'J  l]  will  he  continneil,  concentrating  on  the  efl’ects  of  iniiiioilelh  d  slater  on  I'.n- 
o\(rall  p'l  rurniain  I-  and  rohustiusss  ol  the  estimation  ;nnl  control  system,  1  his  thesis  n^ee,  tin 
Kalman  lilters  in  a  hank  citnflgiirat ion  whi'fe  each  filter  lias  Ix'cn  optiini/ed  for  a  given  coinhinat  ii  iti 
of  the  par;nnelers  of  inti'rest.  One  to  computational  loading  concerns,  the  parame  ters  have  In  ■  ii 
limited  to  two,  ’Flu'  reasoning  hehnul  ihi,s  will  be  develop<,’d  in  ('ha])ter  2,  1  he  robust  ness  analysis 
of  the  controller  provides  insight  into  the  effects  of  using  ;i  n'duced-ordi'r  sitite  filter  mnd.  l, 

/,/  ISitckiirouuil 

'Fhi?  original  fe.-isiliilit y  study  u,sing  a  moving-liank  mnitiph'-model  estimator/conirolli'r  was 
In'  (lent/  [8]  and  ,sub,sequei)l)y  /tresented  by  Maybeck  and  Ib'iitz  [22],  1  In'  Kalman  lilter, 
a  recursive  optimal  filter,  is  the  [iriitiary  building  block  of  the  MM.AE  techniipie,  and  so  it  is  now 
described.  The  develo[)ment  of  the  Kalman  filter  in  thi.s  .section  is  not  intended  to  be  complete, 
Ffir  detailed  discussions  of  the  Kahnan  filter  and  related  topics  .see  [F'']-[20],  1  he  noialion  iiscd  in 
this  '  hesis  frdlows  [IS]  such  that  a  stochastic  (iroci'ss  is  denoted  hy  x  while  a  delerminisi  ic  pr()i'e,ss 
Is  denoted  liy  X, 

/,/,/  Kalman  I'llU.-  The  Kalman  filter  i.s  the  foundation  of  the  iM .\1AK/M M .\C  control 
dosigp  To  use  the  Kalman  filter  approach  to  state  estimation,  viewed  here  as  a  siiecific  form  of 
Hayesian  estimation  [19:129-136],  the  system  of  intere.st  is  a,s.snmed  to  he  adi'qutitely  defined  hy  a 
linear  stocha.stic  state  model  [18:1(53-170],  The  ecpiation  representing  this  sysii'iu  is  given  hy: 


x(t)  =  F(t,)x(/)  -F  B(/)u(/)  -F  G{t)w{t) 


(1.1) 


wi"  ri'  x(  ■ )  r-'i'roS'-n!  s  an  /t-vrclor  si  al<-  process,  ii(  ■)  is  an  /-vectfir  deter  mini. --I  ir  coni  rol  in  pul  ,  F'(  )  is 
an  »d)y-;)  system  dynatiiies  matrix  (which  coiittiiiis  tlu'  parameters  of  interest  m  this  ell'ort ).  B(  )  is 
an  /el>y-r  detertninist  ic  input  mat  rix.  and  G(  )  is  an  n-hy-.s- noi;  e  input  tiiatrix.  I'lie  charact erist  ics 
of  tile  while  (danssiati  noisi'  w  is  of  particular  importance. 

I'he  statistics  of  w(  )  ar('  given  by; 


/■;{w(/)}  =  0  (l.-J) 

/;(w(/)w(/')'^  )  =  -  >')  (I  -f) 

where  Q(-)  is  ati  s-hy-s  matrix  that,  is  s  nniK'tric  and  positive  seniidefinite,  ^(/)  is  the  Dirac  d(  Ita 
fiinction.  Q  n'pri'senis  the  dytiatnics  tioise  strength  l>eing  |>nl  into  the  systi'in.  The  sehTlion  of  (J 
is  oftett  left  to  the  engineering  itit  nit  ion  of  the  desigin'r  [■1;U)0]  tiinl  therefore  will  he  a  cliallenge  in 
tin'  estiniator/controller  design. 

Another  ittiportatit  statistical  property  is  related  to  the  state  process  vi'clor  x.  The  itiitial 
condition  of  x  and  the  certainty  with  which  this  value  is  knowui  is  rerpiired,  I  lie  initial  condition 
mean  value  of  the  state  is  given  by: 

/t{x(to)}  =  xo  (l.l) 

'I'lie  right  harnl  side  of  El(|uatioti  (1.1)  is  the  mean  of  x  at  the  itiitial  time  /q.  Tin'  degree  to  which  the 
initial  condition  is  known  is  given  l,y  ilie  covariance  whicli  is  denoted  hy  Po  and  given  in  Fapiation 
(1.5), 

£'{[x(<n)  -  Xo][x(<o)  -  xo]^}  =  P()  (1-5) 

lo  implement  tite  Kalman  filters,  tlie  discrete  forms  of  tlie  above  eqiititions  are  used.  Kipial  ion 
(1.1)  contains  the  matrices  F,  B,  and  G.  This  same  system  representation  may  lie  expresseil  in  a 
discrete  time  formulation  [18:170-172],  The  state  transition  matrix,  is  derived  using  tlie  system 
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wyiMinic.s  rix  F: 

=  <5(/  -/o)  =  (1  'i) 

it'  F  !.'■  a  coaslam  nialrix.  ^  ilafiiics  how  ilie  stal('  vcrtor  changes  ovar  tin'  s|Mcirn'il  tiiiii'  |><  ikmI. 
A'hlitic. Daily,  the  (Ictarminist  ic  input  matrix  B  may  l)c  I'xprcsscd  in  a  <lisi-iiii-i  imi'  form  ii.sing  <I>. 

Bj{t,)=  $(/,  +  ! .  r)B(r)</r  (1.7) 

(ii.'^rrc'ta-timo  wliite  Gaii.s.'iian  .sy.sti'm  dynamif.s  V('ftor  vvy^  ha.';  a  r()\arianci'  Q,;  that  is  a 

function  of  the  stati'  transi'ioti  matrix,  the  iiois;'  input  m.'itrix,  ;iikI  the  .stn'iiglh  of  the  coni  iiiuoiis- 
time  w; 

Qj(/.)  =  J  r)G(r)Q(r)G^'(r)$'^'(/,  +  ,,r),/-  (l.s) 

III  ti'rms  of  those  Jiscrete-tim<'  malriee.s,  e(|uatioii  (1.1)  is  rewritten  ti.s: 

X(M  =  $(/;,!, _i)x(/,_i)  +  B,/(f,-i)u(/,:_i)-f  Wj(/,_|)  (1.!.)) 

ainl  tlie  time  proi'jagatioii  ec|uatioiis  for  the  conditional  moan  state  esiimati'  of  x  ;uid  the  coc.iriance 
are  given  by  [18: 174- 17.5]: 

x(t”)  =  $(/,-,  t._i  )x(<+_, )  -p  Bd(l._i)u(/,_i)  (1.10) 

P(ir)  =  $(li,/._,)P(<+  ,)$^(t.-,<,_i)  +  Q,i(t._i)  (1.11) 

'[  he.se  equations  give  the  relationship  between  two  consecutive  state  estimates  aiicl  covariaiice.s  at  I  he 
times  t~  and  Note  tliat,  in  order  to  carry  out  tlie  above  calculations,  the  ('stimale  x(t,^_,  )  and 

the  covariance  P(tf_i)  must  be  known.  For  these  to  Ire  determined,  a  measurement  from  the  system 
must  be  taken  and  the  enlin,'  system  of  c(iuation.s  updated.  'I'iic  “-f”  and  superscripts  indicati' 


1-4 


that  the  varial)lc  ol  ii.t,ertst  is  ('valu<ali'il  just,  afti’r  or  just  a  iiifa.'iH'c  iia’iit  is  iiicoi]ii >ra( id, 

ri'siH'ct.ively. 

d  ill-  (liscrolc-l  iiiK'  (or  samph'-iiata)  iiirasiirriiiriil  z  is  (iz-tortiiiiK'd  l>y  a  im  asiin  iiK’nl  jiindt  l: 

z(/,)  =  H(l,)x(t,)  +  y(/,)  (lid) 

I  111'  variahh'  y  is  ;:io  iiusisiirciiii'iit  iiuisc  accouiiliiig  for  tin'  iiiir<‘ii  aiiity  uiili  uliii'li  llii'  ii!'a.-.iiii - 
iiii'iils  art'  oht,<'inii  I.  '1  ho  iiioasiiroiiu'iil  iioiso  is  gonorally  coiisitlorod  to  ho  a  (laiissiaii  disria  l i  -i  i n ir 
proioss  with  zoio  moan  and  a  covariaiicc'  of  R.  1  lie  sysit'iii  modol  dovolniiod  assniin.'s  lhal  (ho 
system  dynamics  noise  w  and  the  measurement  noise  v  are  indi'pt'iuh'iit .  In  addition,  the  noisos 

arc  independent  of  tlie  initial  condition  x(tu).  Once  z  has  lioen  nK'asiiri'd,  tiu'  ontiia'  systt.'in  of 

equations  arc  updated  using  the  follc'.ving  three  equations  (in  addition  to  the  I'lpiat ions  already 
given  in  (1.10)  and  (1.11)  for  the  filter’s  propagation  cycle): 

K(U)  =  P«-)H^«.)[H(/i)P«-)H'''(<.)  +  R-('.)]"‘  (l.llt) 

x(<+)  =  x(t-)  +  K(<,)  [xiU)  -  H(t,)x(/-)]  ( l.M) 

r{tt)^r(l-)-K(l,mt,)P{t-)  (l.lh) 

d'he  bracketed  term  in  liquation  (1.14)  is  called  the  residual  (or  in  some  litoratnro,  innovation.i 
[4,  14])  and  denofetl  by  r(/,  ).  d'he  residual  of  the  Kalman  filter  indicates  how  much  error  correction 
is  required  in  the  filter,  since  it  is  the  difference  between  the  most  recent  measurement  and  the 
liest  prediction  of  that  mea-surernent  based  on  e.stimates  prior  to  that  sample  time,  d'hc  residual  is 
weighted  by  K  and  added  to  the  previous  estimate  of  the  state  to  arrive  at  the  lu'w  estimate  of  the 
state. 

The  discrete-time  Kalman  filter  stochastic  difference  equation  and  propagation  equations  have 
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Imcu  jiresoiitcii  a^-  a  it'vii'w.  Several  important,  nmnbers  in  llie  .sealar  ease,  or  malriees  in  llic  vm  tor 


ea.so.  are  iiivrrlveil  m  these  ecpitU ions.  'I'liey  are  rleriveil  from  tin'  modi'l  of  the  s_\>.lom  ih  vc  lopi’d  hy 


till'  designer. 


I.I.J  System  Model  Tlie  integration  of  control  systi'iii  design  and  strnctnrid  eonsidi-rat ions 
during  tlx.'  design  pha.se  of  a  system  is  becoming  l)otli  commonplace  and  necessary  [iha:  ITiiS],  While 
this  research  effort  doi'S  not  develop  the  actual  system  model,  it  is  important  to  nnderstand  some 
basic  concepts  that  underlie  the  morhds.  'The  large  space  st met nre  considered  is  modelled  by  hii' ar 
ei  piai  ions, 

bhe  thesis  resetirch  performed  over  the  last  live  yetirs  has  involveil  dilferent  types  ol  investi- 
gation.s.  Ilenlz  /lerformed  ;\  feasibility  study  on  using  a  moving-bank  MM.\b/.\l.M.\( '  approach  to 
solve  the  control  prohlem.  I'ilios  continued  to  look  at  the  problem  of  apjilying  this  same  techidipie 
to  ;i  more  realistic  spacecraft  modt'l:  however,  onlint'  adaptation  was  iml  really  rei|nired  for  that 
parlicnlar  siiacecraft ,  Kariiick  was  given  the  ba.sic  two-bay  truss  mmlel,  .uldeil  niodilicai  ions,  and 
used  the  new  model  to  pi'rform  his  researclt.  La-shlee  ii.scd  the  same  moch.'l  to  perform  his  rr'si-arch. 
V'aii  Der  Werken  used  the  same  filter  model  but  included  higher  order  mode.s  in  the  truth  model, 
2-1  states  versus  (3  states. 

One  of  the  iirimary  concerns  in  control  design  is  disturbance  rejection.  'I  Ik'  ahilitv  to  with¬ 
stand  outside  interference  from  natural  and  man-made  phenomena  as  well  as  disinrhances  from  iIk- 
systi'in  itself  is  essential. 

Ill  large’ s|)ace  structures  the  disinrhances  are  tlic  result  of  slewing/point  ing  maneuvers, 
thermal  transients,  and  mechanical  machinery  such  as  coolers,  generators,  etc.  Control 
of  the  dynatnic  resjron.se  is  e.s.seiitial  for  maintaining  the  ride  quality  and  pi'rforniaiice 
reipnremeiiis,  a,s  well  as  for  the  safety  of  tile  structure  [25;  I7t38]. 

The  parameters  that  affect  the  performance  of  the  system  may  hi’  categori/.ed  as  (I)  the  phvsicd 
rom|iosition  of  the  structure,  (2)  the  sources  of  disturbances,  and  (5)  the  control  system  itself. 
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Figur<‘  1.1.  Rotating  Two- Bay  Tru.ss  Model 

The  [ihy.sical  compo.sition  of  the  .structure  is  shown  in  Figure  1.1.  The  structure  has  co-located 
acceleronieter.s  and  tlirusters,  whicli  sitnplifies  some  of  the  calculations.  Tl'e  huh  has  gyroscop.i's  tcj 
deierniiiie  angular  di.splacenient  and  velocity  and  a  co-located  inertia  wheel  as  an  actu.ator.  The 
.sen.sor  and  actuator  locations  and  functions  will  be  explained  in  detail  in  Chapter  3. 

Karnick  presents  a  very  thorough  description  of  the  components  of  the  two-bay  truss  mothd 
th.it  represents  an  appendage  of  a  larger  space  structure  (such  as  a  manned  space  station)  [1  l;39-oS]. 
The  aluminum  rods  that  make  up  the  two-bay  truss  are  assumed  to  be  of  a  certain  cross-sectional 
area  and  elasticity.  The  mass  and  stiffness  matrices  that  describe  the  model  were  calculated  u.^ing 
finite  element  analysis  [11,  26].  Chapter  3  presents  the  results  of  the  finite  element  analvsis  .and 
the  development  of  the  matrices  that  describe  the  system.  The  matricivs  are  cotisidered  to  he  the 
nominal  values  for  design  purposes.  During  this  research  the  paratneters  that  arc  considered  to 
vary,  namely  the  dainping  factor  and  natural  frequency  of  the  structure,  are  the  result  of  ch.ingcs 
in  the  mass  and  stiffness  matrices.  The  trtith  model  is  developed  by  using  twelve  modi  .s,  the  rigid 
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hotly  aiitl  eli'vt'ii  niotlfs,  wliifli  are  represented  hy  twenty-four  slates. 

Ill  elforls  prior  to  tliost'  that  adilri'ssed  the  rotating  two-bay  truss,  the  .\1.M.\M  and  M.M.\(' 
alyt  uit  huts  were  used  to  coiitrol  models  of  various  complexities.  Ileiilz  used  a  simph'  twtj-stale 
nil  del  of  the  siniciiire  [8:1(,!].  'I’he  control  ratio  of  the  structure  dynamics  is 


(■(s) 


R(s)  S-  -t-  ■2C-t-'nS  4- 

or  in  .state  variabh-  form  (standard  coul  rollahle  form)  tlie  structure  tlynamics  are  ih'scrihet]  by: 


1 

X(t)  =: 

1 

0  1 

x(f)  + 

0 

u{l) 


<:(l)  = 


1  0 


x(/) 


(i.Ifj 


(  I.IS) 


'Ihe  state  vector  x  is  composed  of  two  components:  X|(<),  which  reiire.sents  a  position  variable,  and 
x-j(0'  'vhich  re|)reseuts  a  velocity  variable. 

Th('  work  pi.'rformed  by  f’ilios  [():3<S-41]  used  a  more  complex  model,  f  ilios  (Wtiluated  the 
.M.MAE  ami  M.\I.\('  algorithms  against  the  Draper  Laboratory/ Rocket  I'ropulsion  Laboratory  ('ou- 
liguration  model,  which  consists  of  four  cantilevered  appendages  attached  to  a  central  IiuIj.  1  his 
model  more  closely  represents  a  satellite  with  four  “whip”  antennas  than  the  structure  developed 
by  Karnick.  A  six  state  model  was  u.scd  by  Filios. 

The  work  performed  by  Karnick  [11]  used  a  two-bay  truss  model  attached  to  a  hvtb.  Karnick 
usi;d  a  six-state  model  accounting  for  the  rigid  body  mode  and  two  beiuling  modes  of  t  in'  si  riicl  iirr: 
a  rcduceil  order  modi’l  derived  from  an  analysis  using  the  method  of  siugul.ar  perl  iirbal  inns  [hi]. 
Lashlee's  resetneh  [bo]  used  the  same  model.  Doth  tin'  truth  and  liltt-r  modc'ls  used  six  states  for 
Karnick’s  and  Lashlee's  research,  and  thus  lliere  were  no  immodelled  slate  effects,  Reiiig  concerned 
afiout  the  effect  of  iinmcodelled  states  (a  robustness  concern),  Van  Der  W’erken  ii.si'd  a  'it-slale  truth 
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luoili'l  and  a  G-slati'  rilt('r  uu>d)'l  [21].  I  ho  full  dovelopiiH'iit  of  llu'  triitli  mod,  I  is  picsmli  d  in 
(  liaplt'r  2. 

dill'  nmul'i'r  of  stall's  doscriliin;^  a  model  used  in  control  analysis  cannot  lie  ]>redi'li'nirmi  d 
soli'ly  on  the  basis  of  mat  hemal  ical  insight  or  engineering  intuition.  I  he  computational  loading 
rcipiirements  of  a  specilic  implementation  may  preclude  ihe  usi'  of  more  than  a  hav  sl.-iti-s,  1  liis 
limitation  m;iy  lia  \  ,■  im  pi  ical  ions  on  tin-  rohusi  ness  of  the  control  system  and  l  li,'  concein  of  i  h  is 
t  liesis  elfort .  liy  increasing  t  he  numlier  of  St  ates  \ised  in  the  truth  model,  N'aii  l)er  W’erken  ol  isei\ ,  d 
changes  in  the  system  [lerformance  as  com|)are<l  lolh.'il  portrayed  liy  l.ashlee.  I'hese  ol cser\ .at  ions 
indicate  that  the  (liter  model  st.'ite  size  may  require  modilical ions  for  proper  performance  from  tin 

MMAE/MMAC'  algorithms, 

I.I.  J  .\[ultipl(  MoiUt  Aildptirt  [■.'.\liiiialioti  -  d/.l/.l /■.'  I  he  use  of  full  slati'  feedback  is  I  In' 
nor  lin'd  init  ial  appri  i.icli  for  opt  im;d  coni  rol.  f'lill  slate  feedback  proviih's  tin.'  coiii  red  s>steiii  si  abilii  >' 
ainl  robustness  toexti'rmd  and  inlt'rnal  disturbances,  llc^wi'ver,  perf'ct  .icct'.ss  to  all  slates  is  seldom 
attainable  and  .so  a  .st;ite  e.stimator  is  u.sed  to  provide  best  esiimati's  of  plant  stall's  to  the  controller 
that  wa.s  designed  as  tliough  full  state  access  were  actually  available,  d  he  controller  for  this  research 
is  an  bCJC'  controller. 

I  sing  LQ(.l  control  assumes  that,  the  system  model  is  Linear,  there  is  a  (Jiiadralic  .o.^l 
cril  erion  assoc  iateil  with  the  coni  rol  problem,  and  the  iioisi's  u.sed  in  l  he  dynamics  and  measnreii  lenl 
equations  are  taken  to  be  Gau.ssian,  “d’lie  I^QG  optimal  controller  li;is  the  ((rhinili/  (qiiiidlnK  i 
pjoperty'l'IO'.n].  d'he  certainty  equivalence  property  is  a  special  case  of  Ihe  .separation  pro|ierty. 
d  he  l.QG  controller  is  eriuivalent  to  the  optimal  deterministic  controller  "but  with  the  slate  replaced 
«iih  the  conditional  mean  of  the  state  given  theob.served  mea-surements..."  [20:17].  Figure  l.2.shou',s 
the  flow  of  tlie  sigmds  from  the  system  to  the  estimator,  d'he  I'stimatc.ir  provides  the  current  state 
estimation  to  the  conlroller,  which  in  turn  provides  a  control  signal  to  the  systi'in. 

Till'  usual  [iroblems  encountered  in  a  real  world  system,  ;is  is  the  rase  with  the  (h'xible  space 
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External 
Di s  t  urbcnces 


Figure  1.2.  Diagram  of  Feedback  Control  System  [20] 


structure,  arc;  that  the  states  themselves  may  not  be  directly  metusuraldo  or  the  .system  par;uii''teis 
may  vary  from  the  original  design  tvs  discussed  earlier.  The  approach  th.ai  M.M.\F  uses  is  to  de.^mi; 
each  Kalman  filti'r  for  a  disc  rete  point  of  interest  in  thi‘  param>;ter  space',  '['he,'  rt.'al  paitimc'tc  r~ 
may  airtii.illy  be  contimioiis  in  th<.'  region  of  interest.  The  nutttber  of  valiic.’.s  that  a  par.ametc  r 
could  attain  would  theic.'fore  be  inlinite  and  the  resulting  comliinations  of  p.-irtinic.'ters  wciiild  al>.i 
be  infinite.  Since  the  use  of  iufiuiti'  (or  even  very  large)  compntc'rs  is  not  rc'.ali.'iic  or  perac;!  ic.-ii 
for  this  research,  the  parameter  space  is  viewed  as  a  discretized  version  of  this  continuum.  Fi-r 
example,  if  the  system  had  tw'o  parameters  of  interest  that  could  take  on  two  di.screte  values  each, 
then  the  parameter  space  would  be  comprised  of  four  discrete  points.  If,  however,  the  parameters 
could  take  on  100  values  each,  the  parameter  space  would  be  composed  of  10,000  discrete  points, 
each  recpiiring  a  Kalman  filter  and  eventually  an  LQG  controller,  which  would  be  an  unbearable 
cotnputational  load  for  any  real-time  control  system.  The  work  performed  in  previous  research 
[11,  15,  2d]  uses  a  two- parameter,  lO-value-per-parameter  approximation  of  the  actual  continuous 
parameter  space.  The  parameter  space  employs  100  points.  Figure  1.3  shows  the  block  diagram  of 
the  MM.VE  algorithm.  A  full-scale  MMAE  algorithm  would  require  a  bank  of  100  Kalnuin  filters 


for  thi.s  applicatioti 


1  li.'  .i.lai.i  i  \  I  -  '  si  I  iii.il  ',1  m  ( x'CMrs  line  U)  I  Ik'  itiiiIi  iplK'.il  k  iIi  i  iI’  I  hi-  pri  ihalu  1 1 1  k  -  :;f  in  ral  ■  >  I  I ,  isi  ^ 

'll  till'  I'-suliials  f'.i'ui  I'acli  Kaliaaii  (illcr.  The  stalo  l•s(i|llat<'  of  ^•al•|l  filur  is  iiiiill  i,)li'  I  liv  lli<' 

I'l'lialiility  lliat  Ilia  assiinu'il  vatiK'  af  paraiiK'lf'r  for  tlial  pari  iciifir  filler  is  ila'  lu-sl  paraiin  irr 
\  iliia  to  Use.  11a'  pruhalalily  Iroiji  I'acli  Kalman  filter  is  a  fmiclKin  of  llu-  <a>iiilii  i' aial  pn  ilial  i  li  l  \ 
'1'  ' i-ir  i''"  I'l'lal I'll  I ly 

■/^Zl  ( ,  l|il.Zi  !,-i  |(^''  I  1  1  ) 

^ - :  '' 

- 1  -^z(  I,  i|a.Z(  ( I  ‘'j  •  -  1  i-i  ) 

'■'li'i''  K  is  ill''  iiiiiniii'i'  "I'  lilti'i's.  fill'  lirsi  Miiiiii'rali  ir  li-rm  repri'si-iils  ilia  pa  .lial'ilii  \  .|''iisii\  :  y 

ill'  ''iirr  111  iiiaasiii'i'mi'iii  hasail  oa.  tin-  assi.iiiail  a*,  paraiiialar  ami  nivan  ilia  pravii  in,-  i  iiaasn  r. 'i  i  a  nl 

I  i  1  -  ’  ’  'I','  I  li  I'l  'll  i;li  Mill''  _  1  .  Z'  _  '  ■  I  lia  samm!  miiiii'i'ali  ir  lariii  is  i  lia  pr  a  ni'  /',  \  alia ',  ami  I  liii,~ 

I  1  I 1  IS  an  III  I' ii  I'll  I',  'a I  a  III.  f  im  pii ihaliilily  ilansi.y  liim  i  m iii  is  f  .rmi' i  1  .y . 

{■}  =  {-~rUi,)A;'(t,)n{i,}]  (i.'jD) 

will  I'a  A,(f,)  is  a  liiiii'iion  ofllia  maasiiramanl  matrix  Hn-ff,),  the  cot  aria -laa  malrix  Pa(/”).  ami 
'll'  ii"i>a  aovanama  malrix  Oa-IM.  namely  )H(  (/, )  +  HaP, )]  ■  Ilia  'f- 

li'  iiiiiiai  "I-  IS  ilia  ,->111 II  I  if  all  numeral  or  terms  for  j  =  I .  .  .  . ,  so  lliat.  I  lie  simi  <  'I  all  /n-  is  iniit  \ . 
.\ii''llii  r  w.'iy  I'l  sialiii'j,  tins  is  that  the  romlilional  prohaluhty  is  the  prohahilil  v  that  iha 

'h'  I'l'la  raml'iin  pimass,  a,  eipials  samii’  sperilic  value,  aj.,  eiven  that  ihi-  maasnrameiil  Zl/,)  has 

lak'ii  a  spriilji-  I'lalizat  ion  Z,  [  1  h:  1  .'iO- Kl  1],  I  he  weighii'd  I'slimates  are  then  sinnniad  to  provi'la 
'll'  hi  •.!  pi' ihal'ihsi  i'  .'illy  W'iithleil  e.slimale  of  the  rnrrent  stale,  '1  he  rnnipiil  ail  pKili.ihihl  a  s  ha\a 
1 1,"  '  li.'ira'  lariat  !'■  ihal 

K 

=  1  (l.-ifl) 

1  =  1 

I  Ins  property  if  las  not  |el  any  one  (liter  rout  rol  I  he  i-sl  imal  ion  piroce.ss  unless  one  of  t  he  esi  imat  i  ir' e 
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I  '  I  I  1 '  >1 1 1 ;  \  ^  I  I  .  til,  ;i  1  .-,11  ill,  ,  ,1  ||l•|■s  go  I  O  /I  TO.  I  ||i'  |  lO.s.-,  1 1 1 1 1 1 1  y  1  li  ii  1  I  1 1 1>  r;i>, '  i',  ,1 1 1,  I  ,  i,  t  ii  I 

!i  I'  i  '  '  I;  t  ,  i|iM.  |.  1 ,  1 1  |,\  r(>«',T,,|i  I'lForls  [N,  ‘>'1.  tl.  II,  15,  21.  I  Ij,  iiiii-'l  1m’ 

I  'I '  \ '  111  in>  jTt  tl  till  'lilt  If,-'  I  roll  I  '.’.oiiii;  lo  /iio,  siiiri'  any  siii'li  prohahilily  WoiiM  ln-n  .nif  ]Mri  ii.iii'  iil  K 
lot  K' t|  tiiiio  ziTo  lt\  ill.-  iliralioii  jiiTl'orini-d  liy  (1.19)  <'V<'n  if  rhanging  coiiirit  ii  m.-'  w.  ro  t"  fan't 
ill  ii  |i,'uUtul,ir  I  i.i  I  Ml  I  It  I  IT  Ntilui-  It)  ln'conii'  tlio  hcst  valui'.  Oik’  way  to  |ii'i  V'ni  ilif  |  ii,  ,1  tin  ii  I  n  a 
lit  III  gtting  lo  /,  I't  1  IS  lo  si't  solin'  lowiT  ami  uplnT  liinils  that  iIh'  //*-  may  lakr  on.  Om''  ilii-  i> 
,Mt  ■  'iii|  tli^li,  tl ,  a  \  i.i  I  t)i’  aila|  il  i  vi'  si  at,-  I'st  iiiiali'  is  availaMi-  for  nsi'  Ky  a  <-oiil  rol  s\  st  '  iii  as  sla  tw  n  i  n 

I  1  j.iiiM'  1  2  I  ir  1 .1 . 

Ill  MiillijiU  Miiih  I  .\(liii>tiii  ('oiiliol  -  MMAC  Karnick  il'Vi’l.  ip.'i  I  tlir  l,(.H','  im  ini  |t  i||,|  j.  ,| 
lilt  I  wi  I- 1,  ay  I  niss  motlil  [1  1  '.IT-!.)''],  I  iial  lia.s  Ix-i'ii  used  liy  oi  1i<t  rcsi-aiTli  rirori  s  [  1  2  Ij ,  1  In'  .M  M  .\( ' 

■  d  Z'  ril  Inn  uses  ilia  .\1  M  .\  1'.  .il.yoiii  Inn  lo  di'V('lo|)  I'st  imali-s  of  l  li,-  st  ai  is.  I  .atli  o|  i  lif  si  ,it  a  '  si  1 1  nat  as 
Is  I  in  II  |i  1 1  into  a  ft  .m  i,  ill,  |'  ,  ,|,i  nnr/ail  lor  t lie  |>araini‘l<‘r  valiia  forn  siioinliinj,  l o  ii >  asi  ii nal  a  in | ,i]i 
I  I'-ura  I  I  shows  f lia  sigd.il  llotv.  '!  Ii<-  blocks  ')<tk>I<'<.I  liy  —G’f/Ji-)  ara  o/'iiiiial  fonirollar  gains 
s|,' t  ilitMally  ast al ilislii't  1  for  aacli  discrete  paraineti'r  value  </*..  Sliglil  mist  lining  of  I  lie  Idler  and/or 
fonirollar  paraiiiatars  can  cause  |>oor  performance  of  an  individual  LQ(1  conirollar.  'I  lia  M,\l.\(.' 
ff.tii! roller  [irovides  .some  relief  from  this  |)arameter  mismatch  jiiohlein  .since  it  uses  a  weighted  snm 
of  fonirol  inputs  liased  on  the  hypothesis  conditional  prohahilities  computed  ni  the  basis  of  the 
res), Inals  generatetl  hy  the  Kalman  filters.  'Ihi'  MMA(.!  approach  lo  a  control  solution  .suffers  from 
the  same  coinpiii  at  ional  loading  problem  as  MMAK.  An  alternative  to  fnll-srale  M  M.\K/M  M  .At ' 
IS  the  iisi-  ol  moving-hank  .\!  .\1  A  h/ M  ,M  A( !  to  provide  a,  very  close  aiiproximalion  to  fnll-srale 
M  .\1  .\ K/ M  M  .\(  optimal  control  wit  h  reduced  compnlat ional  leading. 

/. /.T  M  on  lit/- Hinik  M  M  A  I'./M  M  A( '  1  he  tln'ory  In'himl  thi'  moving- 1  la  nk  M  \I  .\  K/M  M  .\( ' 
1'  ilial  llie  eslimale.  and  in  turn  the  propir  control  input  to  I  he  system,  max'  be  a|i|iroxmiated 
l)\  a  rednred  nimiber  'if  fli.scr'l  ized  parameter  values  that  lake  on  the  values  of  approiniai  e  (d\- 
iiamifally  rer|er|ared)  local  ions  used  for  t  he  basis  of  lln-  full  bank,  f  igiire  1,.m  sho'vs  how  ihe  full 

M.'t 


Figure  1.5.  Diagram  of  MM.VE/MM.VC 


^I.^I.■\E/.^I.M.AC  system  would  look  if  all  the  filters  were  used.  Each  of  the  blocks  repre.seiit  a  di.s- 
cretij  parameter  point  used  as  the  basis  for  a  Kalman  filter  and/or  a  control  gain.  Figure  1.6  shou'.s 
how  a  moving  3-by-3  bank  might  appear  surrounding  the  current  (estimated)  parameter  value.  The 
ntoving  bank  is  composed  of  nine  solid  blocks.  The  “discarded  filler”  points  correspond  to  a  3-by-3 
bank  at  an  earlier  time  instant  when  that  set  of  nine  points  most  closely  surrounded  the  estimatc'tl 
parameter  location  at  that  time. 

The  problem  of  how  to  move  the  bank  has  been  discussed  and  evaluated  [8,  22,  6,  11,  15,  2  1]. 
Thi.'re  are  five  ba.sic  techniques  for  deciding  to  move  the  bank  or  perform  a  coar.ser  rcdiscretizatinn 
of  the  active  bank  of  filters  as  shown  in  Figure  1.7.  The  techniques  are: 


•  Residual  monitoring 

•  Parameter  po.sition  estimate  monitoring 


•  F’arameter  position  and  “velocity”  estimate  monitoring 

•  Probability  monitoring 
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figure  1.6.  Movirig-bauk  MMAE/MMAC  -  Fine  Discretization 
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Figure  1.7.  Moving-bank  MMAE/MMAC  -Coarse  Discretization 
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•  Parameter  esllmatiou  error  covariance  monitoring 

R.osicliial  inoiiitoriiig  uses  a  likelihood  quotient  delined  as; 

U.{ii)  =  ii(/,)A;'i-,(/,)  (I 

t  ■  vieoi,inuue  wncther  a  move  ot  the  u«nk  is  require<i,  wiiere  Aj'  l.>  «itt:  ii.Vci'.-vL'  ■’ I  *1  ‘ 

in  liquation  (l.l-'i)  (liie  bracketed  inverse).  If  all  the  Lt  are  above  a  preset  bound,  the  bank  Is 
mo\ed.  Furtlier,  tlie  filter  with  the  smallest  likelihood  quoliiMit  should  be  tlu'  lilii-r  nearest  the 
true  parameter  value,  llesidutil  moiiiloring  is  susceptible  to  singh'  large  sam])les  of  measurement 
noise  and  may  give  false  alarms  [22:92],  i.e.,  the  bank  may  be  moved  evt'ii  I  hough  the'  move  is 
umieccessary. 

Parameter  Position  Estimate  Monitoring  attempts  to  kee])  the  btink  centered  on  the 
current  estimate  of  the  true  parameter.  The  estimated  true  paranu'ter  is  given  by: 

K 

nil,]  =  £’{a  i  Z(ti)}  =  ]^a|.  pki!,)  (l■2.■i) 

*•=1 

When  the  "distance”  between  the  center  of  the  bank  and  the  estimated  parameter  loctitioii  becomi's 
larger  than  a  preset  value,  the  filter  bank  is  moved.  Since  the  calculation  dcpetids  on  the  time  history 
of  measurements,  this  technique  is  less  susceptible  to  false  alarms  thati  the  residual  monitoring 
method  [22:29]. 

Parameter  Position  and  “Velocity”  Estimate  Monitoring  is  an  e.xtension  of  the  |>re- 
vious  method.  By  tracking  the  “velocity”  of  the  parameter  estimate  through  the  most  recent 
parameter  position  ('st  irnates,  the  ne.xt  position  of  the  parameter  may  be  predicted.  If  the  dislanci' 
between  the  predicted  location  and  the  current  center  of  the  filter  bank  e.xceed  a  preset  limit,  the 
bank  is  moved. 

Pi  ■cihfihility  Monitoring  via  t)ie  eqiiatioii-s  already  pr<'.sent<'<l  also  jirovides  an  indication  ol 
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till'  most  likely  location  of  the  parameter.  The  computed  prohahilit ies  are  compared  to  a  ineset 
threshold  or  simply  the  largest  prohahility  to  determine  if  the  ci'iiler  of  the  hlter  hank  .should  he 
moved  to  the  filter  that  is  |irodiicing  the  probability  exceeding  the  threshold. 

Paraniotor  Estimation  Error  Covariance  Monitoring  provides  a  means  for  deterniimng 
whether  the'  liank  should  be  contracted  from  a  coarser  discretization,  as  shown  in  Figure  1.7.  I'he 
parameter  estimation  error  covariance,  given  by  Fijnation  (1.2-1),  i.s  monitored.  I'he  jiarameier 
eo\ari:mce  can  be  readily  cotninited  as: 

P.d'M  =  /■;{[a-a(/,)][a-a(/,)j'  I  Z(/,)  =  Z,} 

K 

=  ^[a  -  a(/.  )][a  -  a(/,  )]'^'  ■  w-f/, )  ( 1 .2  1; 

*  =  i 

within  the  MMAE  or  .\IM.\C'  algorithm.  W'hen  the  norm  of  the  vector  P„  falls  lielow  some  inesei 
threshold,  the  batik  may  then  be  contracted.  Monitoring  the  parameter  estimation  ei  roi  covariance 
is  tiot  effective  for  expattsioti  decisioti  making,  however,  since  Eqnatioti  ( 1.2-1)  depends  oti  the  ctirretit 
choice  of  aj.  atid  is  "artificially  bounded  from  above  by  the  current  size  of  the  baid-;"  [22J. 

W'heti  the  |iarameter  undergoes  a  jump  cliange  or  i.s  changing  rapidly,  the  clo.soly  spaced 
filter  bank  tnay  tiot  be  able  to  track  the  parameter  adeipiatcly.  Hy  exiianding  the  filter  batik,  the 
[laramcter  is  reacquired  and  a  decis’ion  to  contract  the  bank  around  the  iiew  parameter  location  can 
be  made.  Rcturnitig  to  Eqtiatioti  (1.22),  it  is  seen  that  this  eijiiation  is  depi'inh'ut  on  the  residuals, 
r(..  If  all  the  likelihood  ipiotietits  ari'  large  and  close  in  value,  then  a  good  estimate  of  the  parametc-r 
location  in  the  paratneter  s[)ace  is  not  jiossible.  Since  the  estimator  does  not  "know"  where  the 
[larameter  is  located  with  any  certainty,  it  is  more  approjiriate  to  exp.atnl  tin  bank  rather  thtin 
move  it  [22], 

licntz  found  that  probability  monitoring  is  the  method  of  choice  [22:9.a].  bashlee  used  the 
same  technique  atid  observed  the  effects  of  varying  the  |)arameter  sjiace  discrel  izal  ion  oti  the  ability 
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of  the  filter  to  track  the  true  parameter  location  [15],  Botli  of  tlu'se  efforts  i)ri’seut  an  extensive 
discussioti  of  the  decisioti  logic. 

'The  paraineti  r  est  imate  of  the  moving  bank  MM  A  b/M  M  A<f  algor i t  lim  should  lie  as  neai  ly  as 
good  as  the  full  hank  I'stimate,  jtrovided  that  the  va.st  majority  of  tin'  full-scale  .MMAK  tilgoriihm 
partimeter  probability  weight  is  contained  witiiin  the  moving  bank,  btit  problems  have  been  sei u 
in  five  years  of  thesis  research  [8,  22,  b,  11,  15,  2-1].  The  tuning  of  the  noise  matrices,  Q  and  R. 
the'  discretization  strategy  used  to  generate  the  discrete  paratneter  points,  ;ts  well  as  the  tuning  of 
the  other  parameters  in  the  filter  algorithm,  apparently  affect  the  system  rather  snhstatil  iaily.  \’;m 
I')'!-  W'erkencs  thesis  shows  that  the  number  of  states  used  in  the  state  s])ace  represent  at  ic  m  ol'  the 
truth  model  may  be  a  very  important  cause  of  poor  tracking  ])<,'rformance  of  the  M  .\1  .A !',/ M  .\1  ' 

used  to  control  thc'  spttce  structuri'  motiel,  and  so  this  robustne.ss  i.ssue  w.irrants  furl  her  siud_\, 

J.2  Frobkid 

The  computat iotial  burden  of  a  fixed-bank  MMAE/MMAC  implementation  is  higher  than 
a  moving-bank  eslimator/coutroller,  and  so  the  rnoving-bank  algorithms  are  preferable'  for  imph'- 
menfation,  Thc  problem  remaiits  to  determine  the  robnstnass  of  tlic  estimator  aiul/or  controller 
to  unmodellecl  effects,  wliether  it  be  a  moving-bank  algorithm  or  a  fixed-bank  algorithm  tising  a 
necessarily  coarser  parameter  discretization  to  maintain  computational  feasibility.  I'licrefore,  the 
mismatch  of  the  2'1-state  truth  model  and  the  6-state  filter  model  is  investigated,  'f  he  selected  truth 
model  is  a,s.surned  to  model  the  real  world  adequately.  This  thesis  re.searrh  will  provide  insight  into 
the  effects  of  the  mismatch  and  determine  the  number  of  filU'r  states  required  to  track  llu-  states 
and  control  the  space  structure  adeipiately. 
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I.  I 

The  researcli  pcrforined  investigates  the  control  of  a  large  space  struct  nre  r('pres('nte(l  l)>' 
a  rotating  two-l)ay  truss  as  developed  by  Karnick  [11].  The  Iru.ss  is  allowed  lo  move  only  in 
the  \-y  plane  of  Figure  1.1.  Restrict ing  the  motion  lo  this  i>lane  allows  for  the  investigation  of  the 
.\IM.\F/.\1,\1AC  algorithm  without  overly  complicating  the  model  witli  motion  out  of  this  plane,  i.e, 
without  coupling  ijf  x-y  plane,  x-z  plane,  and  torsional  motions,  Masses  attachcsl  to  the'  structure 
were  added  by  Karnick  to  allow  changes  in  the  mass  matrices  lo  obtain  a  mori'  realistic  sp.ice 
structure  model.  This  situation  represents  the  effects  of  depletion  of  fuel  tanks  on  the  structurt'. 
'Fhe  model  developed  by  Karnick  is  described  in  terms  of  ma.ss  and  stiffness  matrices  ohtaiuetl  fri.iiu 
a  finite  element  analysis  [26]. 

The  paraiiK'lers  of  interest  for  this  thesis  are  the  mass  and  sliffne.ss  malrici's  that  a|'pear  in 
the  system  dynamics  equations  through  modal  analysis  [25,  11],  Tuning  of  the  Q  matrix  describing 
the  uiicertainties  ti.ssociated  u'ith  these  dynamics  model  matrices  is  continued  from  the  elforts 
of  Karnick  [11],  btishlee  [15],  and  N'an  Dcr  Werken  [24].  .Additional  tuning  ellecls  due'  to  the 
measurement  noise  covariance  matrix  R  are  also  considered,  d'lu'  mi'asnrement  noise  cmuriance 
matrix  was  considered  to  be  diagonal  in  previous  theses  [15,  24:94,77].  .A  diagonal  matrix  indicates 
that  the  measurement  noi.ses  are  independent,  i.e.  the  measurement  of  the  velocity  at  node  1  is 
not  affecti'd  by  the  measurement  of  the  velocity  at  node  2.  R  will  be  modified  to  include  non-zero 
off-diagonal  values  that  are  cho.sen  in  a  physically  motivated  manner. 
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l.-i  Approach 

Previous  thesis  cITorts  concentrated  on  observing  the  estimation  and  control  of  the  fillt'r- 
modeled  states.  The  filter-model  states  used  in  the  previous  research  effort  are  di'fiiied  by: 


Rigid  body  mode  position 

X-j 

First  bending  mode  position 

Aa 

Second  bending  mode  position 

■'•■1 

Rigid  body  mode  velocity 

XTj 

First  bending  mo<ie  velocil  \ 

J'C, 

Second  bending  mode  velocity 

The  truth  model  vi'ctor  is  correspondingly  assumed  to  be  compostal  of  the  rigid  body  and  first 
eleven  bending  mode  position  states,  followed  by  the  corresponding  twelve  velocity  states: 


■ 

Rigid  body  mode  jiosilion 

First  bending  mode  position 

J'12 

Eleventh  bending  mode  position 

■fia 

Rigid  body  mode  velocity 

First  bending  mode  velocity 

Eleventli  bending  mode  velocity 

,‘\n  error  vector,  r;,  was  used  by  Van  Der  Werken  [2-1]  for  evaluating  e.sti;nator  |n'rformance,  given 
by; 

e  =  Tx(-x/  (1.27) 
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lilt'  suliscriiits  iiulicate  wlit'tlier  tlio  variable  is  a  filter  variable,  f.  or  the  triilli  iiiodi-l  variable,  I. 
rill'  matrix  T  is  a  (j-by-21  transformation  matrix  which  moves  the  states  repri'si'iil  ing  tin'  vel<n  il\ 
of  the  rigid  body  mode  and  the  first  two  bending  modes  to  the  fonrili,  liftli,  and  sixth  vecinr 
liositions.  van  Der  W'orken  compared  the  truth  and  filter  slates  (an  estimate  in  the  filter  case) 
dii'  Ctly.  thus  askitig  whether  the  filter  estimated  the  Initli  modi'l  states  a.s.sociated  only  with  the 
first  three  modes.  'I'liis  tlii'sis  effort  will  instead  ask  if  the  reduced  orrler  filter  can  estiiiiati'  (or  if 
tlie  reduced-order  regulator  can  coittrol)  the  true  total  shape  of  tin'  truss  (ilne  to  all  iiiorh's'  effects). 
Furthermore,  it  will  attempt  to  answer  robustness  cpiestions,  i.e.  will  tin'  controller  work  with  a 
ri'diici'd  stale  filter  model,  by  gradually  increasing  the  higher  ordi'r  state  effects  into  the  system. 

f  li '  new  erri.'i'  term  will  lie  defined  :>s  tin’  differenci'  betwien  true  and  esiimaleil  p's-iliMUs 
and  vehicities  at  tin,'  acceli'ronieter  ;ind  thruster  (used  tis  conlrid  eleinenls)  localiiins  on  the  truss, 
and  at  t  In-  intb,  of  Figure  1.1: 

e  =  H(X(  -HjXj  (1  .L'S) 


where 


H/ 


H/ijxj  Oax.-) 
03x3 


c.x'i 


(1.2S)) 


H,  = 


H 


/  i  3  X  1 


AH 


^  3  X 


(l.liO) 


0 


0 


^/2.vx  .1  I* 
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where  H(p  corri'sponds  to  the  effi'ct  of  the  upper  nine  modal  position  states,  and  H,,  to  tin'  ('fferi 
of  the  iip[)er  nine  modtd  velocity  states,  H/  will  provide  a  means  to  observe  the  elfecis  nf  i  he 
nimiodi'lled  states  in  the  operation  of  the  ('stimator/controller.  I'lie  scalar  mull  iplier  A  will  fie 
varied  from  0.  the  case  where  the  I, ruth  model  is  of  the  same  orrh.'r  as  the  filter  model,  tu  1,  the 
case  where  tin.'  full  eflect  of  all  the  states  of  the  truth  modi'l  are  bi'ing  incor|)orated.  liy  gradually 
increasing  A,  the  effects  of  the  unmodelled  states  sliould  become  aiiparent  .  d'he  rerpiirement  for 
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adding  iiicirc  stairs  to  th<‘  iillcr  model  is  dt'lermiiK'd  l)y  liow  large  A  Ix'eonies  for  tin'  coni  rcilled 
systi'in  to  heccnne  nnstahle  or  exhibit  other  unacceptable  behavior.  Afti'r  (plant  if'ying  t  in'  eirects  of 
the  nnmodelh'd  states,  the  H,,,  and  11,1.  matrices  may  be  further  partitioned  to  determitie  which 
sttites  are  needed  for  a  rolnist  cont rolh'r. 

rile  robustness  (.if  the  ostimator/controih'r  is  the  primary  concetti  of  this  thesis,  holiiislness 
is  the  ability  of  the  I'st imator/controller  to  function  properly  in  tin;  face  of  unnio(lelle(l  efl'eets  dm- 
to  reduced  order  lilter  models.  Slowly  varying  parameters  and  jntiip  changes  iti  tin'  |)arameters  of 
interest  will  he  explored  a.s  well  as  stationary  parameters.  If  the  addition  of  more  states  is  indicated, 
the  work  perforttied  by  previous  thesi.s  efforts  will  lie  re|)e;iled  witli  a  ne'v  filler  struct  nri.’  basi.'d  on 
an  appropriately  higher  dimeii.sional  state  vector. 

’file  Qf  :ind  Ry  matrices  will  be  varied  in  siicli  <a  w;iy  ih.ii  it  nsefnl  inning  strategy  can 
b('  csttdilished.  I’la.-vion.s  efforts  concluded  that  the  tuning  of  these  matrices  are  critical  to  the 
robii.stness  of  the  cunt  rolh.'r  [15,  2'1],  'I'he  previous  research  also  concludi'd  that  the  LQ  la.'gnlator 
wi'ighting  matrices  are  also  critical  to  the  performance  of  the  controlh'r.  d’he  prc'vious  ri'si’arch 
n.sed  a  tuning  proci.css  that  may  not  be  practical  in  all  cases.  Qy  was  held  fixed  while-  R,;  was 
varied  until  best  porformaiice  was  obtained.  Then  R/  was  fixc'd  and  Qj  was  chttnged.  Oin'c  lin- 
values  were  fixed  in  this  manner,  the  remainder  of  the  research  was  pcnformed.  \’an  Der  W'erken 
exiierienced  difncnlty  in  the  tuning  process:  Tlj  was  appropriately  e.stablished  a.s  R.,  plus  additional 
terms  to  account  for  the  impact  of  tlte  18  iinmodoHed  states  (this  will  be  discus.scd  in  detail  in  the 
following  chapters),  but  Q,  was  incorrectly  equated  to  Q/  as  Qj  was  varied  to  accomplish  tuning. 
A  way  to  perform  true  tuning  of  Qy  and  Ry  shall  be  incorporated  into  the  software,  and  the  tuning 
philosophy  of  Qj  and  Ry  will  be  researched.  In  addition,  once  the  noise  matrices  have  bi'i’ii  inneil, 
the  process  of  tuning  tlie  bQ  regulator  weighting  matric<'s  will  ;dso  be  I'valuated. 

Investigation  of  tin'  difliciillii's  encountered  by  itrevious  thesis  resi'arch  in  tuning  tin'  wt'ighling 
matrices  will  be  continued.  Once  the  weighting  matrices  liave  been  tuned,  the  iiarameti'rs  will  he 
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all'.nvt'il  to  traversi;  I  ho  [laraiin'tor  space  slowly  to  st\Kly  the  ahility  oft  lie  M  M  A  h/M  M  AC  al"oiil  Inns 
to  inaiiit ain  coiit  rol .  Jump  changes  in  the  paranieli'rs  will  also  he  invest iga I e<l.  '1  he  para i net er  space 
iliscretizal  ion  itself  will  Ix'  evaluated  on  the  basis  of  bashlee's  recoinmetulat  ions  [15]. 

Laslili'e  oliserai'd  that  the  parameter  idcnt ifving  cajiahility  of  the  estimator  was  enhanced  by 
proper  discretization  of  the  |)aramet('r  space  [15:115-117],  'The  choice  ol  parameter  disciet  izat  ion 
levi-ls  is  determined  by  holding  oik'  parameter  constjuit  and  varying  the  other  parameter  in  tin 
Irnlh  model,  from  the  data  gathered  from  a  single  filter  from  the  bank  matched  against  the 
changing  truth  mode],  rms  I'rror  ph.ils  .are  ii.sed  to  determine  the  disrri't  izat  ion  Icvils  that  yiehl 
a  givi'ii  pi'rcetitage  di-gradat ioti  iti  tin'  likelihood  ([uotietit  givini  by  laination  rl.’J'J),  1  hiis.  the 

rlisci'et  izatioti  is  established  such  tlnit  the  elemental  filters  (and  controllers)  iti  the  liank  are  readil>' 
distitiguishable  from  oite  atiother  ott  the  basis  of  their  residual  clniracti't'ist  ics. 

/.)  Suiiiinnrij 

’['he  remaimng  chtipters  of  this  thesis  covi-r  the  following  areas.  Chapter  '2  de\clops  the 
algorithitis  iisi.’d  in  the  .\1M.\E  attd  MMAC  portions  of  this  re.si'arch.  iticlnding  decision  logic  for 
ttiovitig  attd  expanding/cotitractitig  tint  bank  of  (ilters/controlh'rs.  The  modelling  of  t  he  large  ;~pace 
structure,  the  two-bay  truss  structure,  is  developiul  in  Chajiti'r  5.  ’!  he  simulations  peiloriiied  by 
this  thesis  will  be  e.\|))aiiied  in  Chapter  4.  Chapter  -5  presents  the  results  of  the  research,  with 
corresponding  conclusions  and  recotntnendalions  presented  in  f'ha|iler  (5. 


II.  .\lijorillmi  Dt  n  loi)iii(  III 


/  I  n  1 1  mini  I  it<  II 

I  II'  lU'  il  l  I  li  I  ti  i  li  \ '  li  i|  iiiirhl  fur  t  li>'  com  rol  of  ;i  lar';i'  ;,i  riici  iii-.'  ii>iii'j,  M  M  A  I  ’ /  M  M  \( 

li:m  iii'i  I  liaio;.  I  '](i|iri‘i  iaM_v  over  tin-  last  three  thesis  eirorts  ;ierfV)rnii<l  li^  Karnick  [II  l.'i-.is 

I.  ishlei'  '  I  a.  1  It- I'l  I .  aihl  \’aii  l)er  WVrken  [2-l;2.‘t- l.'l],  'I'his  chapter  [)resents  the  (IrViloiiUieiil  "I  I  li 
al'^i  irit  Inns,  with  liith-  iiioilificat  ion  I’roiii  the  previous  efforts.  .\ew  issues  aiul/.ir  issias  rniir.il  t 
this  thesis  research  will  he  stri'ssed. 

J. J  ISiijjtsiiUi  .Mnlliplt  .M  dill  I  Ailiijilii'f  ll.il  till  lit  Inn  .Miinrillnii  l)i  ri  Idjiiiii  n  I 

riie  hasir  i|evelo|inieiit  of  the  M.M.M'i  algorithm  was  pre.senied  in  Chapter  1.  I  he  r'llluwin 
maliritil  is  presented  for  completein.-ss.  .A  thorough  tretilinenl  of  this  suhii''-t  is  pi-iseni ■  I  h 
MayI.eck  [I':,>:l2'.)-l;il.i]. 

'file  sysleiti  under  iioesl  igal  ion  is  considered  to  he  ;  deipiately  ni'  idelh  d  1 1>  lie  li nea r  si  i  -  ! i.is 
lie  'll  If'  r'  'lice  I'i  |||:|  I  I'  III 


x(/,  +  i )  =  <!>(/,  + 1  )  +  B,/(/,)x(/,  )  +  G,,(/,)w,((/, ) 


frem  wliirh  measur'-jiieut.s  tire  taken  hy  .sampling  via  llie  linear  relationship 


z(/.)=  H(/,)x(/,)-|-y(/,)  (2.2 

.\s  slale.l  III  (  htipli  r  1.  w,;  tiiid  Y  tire  a.ssimied  to  he  iiidependenl  of  etich  ot  le  r  tiiid  /er,  i-nieai 
w  hile  ( iau.ssitiii  disr;--le- 1  idle  noi.se  jiroce.s.ses,  ']  lie  .si at  i.sl irs  of  ilil eresi  tire; 


wli.r.’  Q  i(t I )  IS  |H  i>ii  i\i-  stMiiuli'iinitc  aiui  R(/i )  is  posit  iva-  dcfiiiilc  fiM'  all  f , .  I  Is’  mil  ;al  cniiilii  imi  i  if 
ilh  >tali'  In  in'I  11.  lailly  kisawii  with  ahsoluti-  cortaiiily  ami  is  <l<■^.^■rill<■^l  l.y  iii.  aiis  ol  a  (laii»iaii 
raihloiM  \ial.ir  that  is  assuimd  to  l>o  imlcpeiult'nt  of  ami  y.  1  In'  imaii  ami  .a  j\ ai  laiir.'  ol  lli' 
uiilial  slate  ci  .|i.  1 1 1  a  m  an  gu’en  hy; 


/^{x('())}  =  Xi> 

/'.’{[x(^i)  -  Xoj[x(/,i)  -  Xu]'  }  =  Pu  I-  I" 

I’he  aliove  (  iiiiat  ions  prov  ide  a  starling  point  for  oliserving  t  he  system  of  interest  [lt):7il],  1  la 
aljiive  e((iiafions  liav*'  lieen  presented  willi  little  ex|)lanat ion  of  the  actual  variables.  1  he  variables 
of  interi'st  to  this  p(hnt  are  [ba,  2'1]; 


•  x(/,)  =  a-diirn'iisiiaial  stall'  vecti.ir 

•  ‘t(bf  1  .  bj  =  state  iriiisilion  matrix  which  transitions  the  stale  from  I,  to 

•  ii(/,)  =  t'-.limeiisioiial  determinisl  ic  input  vector 

•  n.db)  coi.treil  input  matrix 

•  w,j(t,)  =:  .s-dimeiisional  while  (lanssian  dynamics  driving  noisi'  vector 

•  G,/(/,)  =;  noise  input  matrix 

•  /,(t,)  =  m  dimensional  measurement  vector 

•  H(t,)  =  measurement  matrix 

•  y(/,)  =  n/aliriieiisioiial  while  (I'anssian  measurement  noise  vector 

f  rom  .Maybeck  [ I '.t:  1  20- 1  .'id],  let  a  In'  the  vector  of  uncertain  p-irameters  (nr  the  model  nmh'r 
simiv.  file  uncertain  iiarameters  affect  any  or  all  of  the  matrices  described  .above,  'fhe  luirpose  o| 
the  liayi'sian  estimator  is  to  compute  the  conditional  di'iisily  function: 
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wliii  h  IS  I  Ilf  Cf  If  lit  ioiial  ilf  iisity  of  x  and  a  given  t  he  iiieasiirei'.ient  history  t  lirniii;li  I  he  eii  ri  eiit  1 1 1  le- . 


Z, ,  where  Z,  is  eoiiiposeil  of  part  it  ions  equal  to  the  realizat  ions  of  z(/ j  ) .  z(  /  j ) . zfl,)-  Has' d  on 

tie-  assiiiiifd  tno'l.  I.  the  lirst  dt'iisity  on  the  right  hand  side  of  Idiuatioii  ('-'.7)  is  (laiissiaii.  I'.iih 
lie  an  xi/^)  and  '  (aarianee  lV>r  a  sperifif  value  of  the  parameter  \ei'lor  a.  1  h'  se  \arialdfs 

ar-  fompiited  l\\  a  l\a!maii  lilter  for  that  particular  parameter  vi-ctor  ;i. 

1  he  iiarameier  vector  a  rati  assume  values  in  the  continuous  range  deliiied  hy: 

,4  C  If' 

u  le  I'e  h'l  IS  real  l.iicli'lean  p-<| I ineiisional  space.  The  )iarameier  c.  i  i.ii-  i,ia\  I"  '  iiii'  ' 1 1  am  l  ail 
C'  li'lani,  It  may  he  slowly  varying,  or  it  may  uiidi'rgo  jump  changi's”  [hai'jl].  Ilf  naiim  ol  a 
I"  iiig  coiilimious  would  causi'  a  jtrolilein  comput at iotially  since  the  calculations  to  s"i\e  hipiaiioii 
(■J.7  j  would  require  an  iiiiiiiite  iiumlier  of  sejiarate  Kalman  filters  ami  int<'grat  ions  that  would  make 
online  usage  of  the  Hayosian  estitiialor  prohibitive.  To  allow  online  computation,  the  parameter 
sp'ace  is  discretized.  The  uliccrtaitt  parantcter  vector  ruiglil  lie  reason;il)|y  defined  as  taking  one 

of  the  valties  ill  lilt,'  finite  .set,  {ai.ag . h/v'}'  wht're  each  discrete  valin.'  is  idii.isen  I'or  optimal 

performance  of  tin,’  control  system.  The  discretization  of  the  [lar.ameter  spact'  is  one  of  tin'  subjects 
tif  this  thesis  effort.  Kach  discrete  value  a*-  has  an  individual  system  moihd  associated  with  it. 
requiring  recom[)iil ation  of  Kqiiat ions  (2.1)  and  (2.2).  Before  the  calculation  of  the  comlit irmal 
m'-an  and  covariance  of  the  Tth  system  state  is  presented,  the  liyjiot he.sis  conditional  prohalality 
will  he  defined. 

I  lf  hypothesis  conditional  probability  is  defined  by 

p*.(/,)  =  pco/i{a  =  a*.  I  Z(/,)  =  Zi} 

_  /z(f.)|a.Z((._,  )(*>  I  1 ) 

Cj  =  i  /z(i,)|a,Z(t, _,)(*'  I  Z,_i)p;(t,_i ) 
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wliirli  rsliow.s  tliai  tlie  conditional  prolitilnlit y  of  the  A-th  filter  being  the  niosl  correct  filter  is  the 
|iioilucI  ot  the  prolitihility  ih'iisity  of  the  most  riTent  measurement  (for  ;i  given  discrete  assniiieil 
paiameter  value  and  tin'  oliservial  measiirinnent  liistory)  and  the  oreviuiis  conditional  inofaliilil  \ , 
divided  hy  the  sum  ol  tin'  A  iirohahility  ih'iisity  aiul  conditional  |"■ollal)ilit v  pi>.i|nci>  (ilnis  I'.ji.-ini; 
the  ;n  (/,  )'s  to  sum  to  otn'  for  all  /,). 

1  he  state  ('stimate  produced  hy  the  estimator  based  on  tin-  ;isstim|)t ioti  that  is  the  .•orrect 
par.atneter  vector  is  given  hy  [19;1;U]: 

x(/  +  )  =  £{x(/,)  I  Z(/,)  =  Z,} 

f  ^ 

—  i  ^ -^xii.iia  Ztt.i^^  I *  '*■' 

U=i 

l< 

=  ^  xrdff  )p*-(<,)  (2.11)) 

k=l 

/he  .state  cst'iiiip.tf'  is  tlte  sum  of  all  tite  probahiJi.sticaily  vvi'igJjted  t'stimates  generated  liy  the 
A'  K'alitian  filtt'rs.  where  the  hy|)othesis  conditional  prohahilities  /<(■(/,)  are  the  weighting  factors, 

1  he  cotnlit iotial  covariaina'  ol  the  sttite  is  calculati'd  in  a  similar  lasliion,  again  weighteil  l,v  ilie 
hyi  )(jtliesi.s  cotidit itnial  probabilities,  and  is  given  by  [lOtlifl]: 

P(/+)  =  £{[x(/.)-x(<+)][x(<,)-x(<+)]^|Z(t,)  =  Z.} 

K 

=  I]pt-(/.){Pfc(<+)  +  [ir.(t  +  )-x(/+)][x,(t+)-x(t+)]^')  (2.11) 

k  =  l 

I  Ilf’  roiidit ioiiai  rovarianrc  Pf(^+)  IS  the  “.stall'  error  covariatice”  associatc'd  with  the  Kaltiian  filter 
iiased  on  the  [jarameter  a^..  'I'Jie  conditional  covariance  is  dependent  on  tin'  measnrenieni  history, 
as  shown  in  Faination  (2.11),  atid  therefore  is  not  precompiitablo.  However,  neither  is  it  absolutely 
re.|iilied  for  onlitie  use  of  the  MMAI)  algorithm. 

'I'he  muliiple  model  adaptive  estiviation  algorithm  developed  above  is  an  adaiitive  filter  atid 
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till'  slructur('  was  shown  in  Figiirc  1.3.  The  discretizalion  of  (he  continuous,  pnrrnnclcr  spncr  into 
a  liiniteil  I'.niiilier  of  points  innkes  online  nse  of  this  algorithm  feasible.  A.s  each  niea.^ntcment  z, 
becomes  available  at  time  the  Kalman  filters  generate  residuals  which  ari'  nsetl  to  generate  new 
by,  lothesis  colulil ional  prcbaliilit les.  With  this  information  availabh',  laiZii  I  ‘*i.Z,._i) 

needed  ill  Kcpiatioii  (2.9)  may  be  computed  for  all  A' as: 


4t'.) 

{  } 


- - r  ‘'^1’  i  l 

(2t)-  I  A,(/,) 

{-^rf(/,)A,T’(/,)ri.(b)) 


where  the  Atli  Kalman  filter  gi'iierate.s  a.s 


(2.12) 


A*.(/.)  =  H*.(T)Pit-(/r)H[(/,)  +  n.r.(/,)  (2.13) 


I'ho  A'  imphimented  filters  thus  provide  the  terim  to  evaluate  the  eonditioiia)  probabiJilie.s  rei/uirod 
to  determine  the  state  of  the  system  given  the  current  measurement.  Using  the  weighting  techniijue 
seen  in  Equations  (2.10)  and  (2.11),  the  state  estimate  x(t^)  is  generated. 

The  othi.'r  (.'stimate  that  may  be  desired  or  recpiired  is  that  of  the  paiaiiii'l er  vretor.  I'lie 
conditional  mean  of  the  parameter  vector  at  time/;  is  given  by  [19:132]; 


a(f.)  s  £7{a(<,)(Z(/i)}  =Z. 

=  /  «4|Z((.)(«  I 

•/  —CO 

K 

=  (2  11) 

A:=l 


'1  lif“  covaririncf'  of  tlic  pstirnatecl  parametor  vector  gives  an  indication  of  tlie  precision  of  the  estimate. 


aiivl  it  can  l')0  rstimati’d  via 

K 

£'{[a-a(/.)][a-a(^)]'''  j  Z{/.)=  Z.}  =  ^fa*  -  a(/,  )][a,  -  a(/, )]  V(/, )  (2.15) 

k  =  l 

It  .should  he  iioteil  t.liat  neither  the  raloitl'.'ions  the  ;iarain<'ier  ostimali'  I'-ir  th'-  covariaiu'-  <>! 
the  parameter  estimate  are  recpiircd  for  the  calculation  of  the  state  I'stimate 

For  tlie  case'  of  non-varying  [)aranieters,  the  MMAF  has  been  shown  to  he  o|)timal  and 
converge  [7],  rin'  situation  examined  hy  Hawkes  and  Moore  sliowed  tliat  the  MMAF  convergeil  to 
the  tilter  that  tiiatched  the  actual  parameter  or  was  closest,  to  the  true  [larameter  [7].  'I  he  [u-evious 
theses  did  not  (le\eloj)  thf.'oretical  restilts  for  the  case  of  varying  parameters  [il.  S.  11.  15.  2  1],  One 
of  the  conci.Tiis  presenti.'d  hy  the  previous  theses  is  that  the  (ilt<'r  may  convergi'  to  one  parami  ii-r 
value  and  lock  on  to  that  value.  If  the  real  parameter  remaitis  at  that  location,  tin’  filter  would 
|)erform  correctly.  However,  if  the  parameter  is  likely  to  vary,  the  filter  would  he  of  little  use  in 
providing  a  state  or  parameter  estitnate.  Dasgupta  and  Westphal  showed  that  the  algorithm  may 
cotiverge  to  the  wrong  parameter  value  for  the  case  of  unknown  biases  in  the  measurement  process 
[5:014-615]. 

To  [jrevctit  the  problems  di.scus.scd  above,  the  addition  of  driiamic.s  p.st.’udonoi.se  t(.)  the 
tiKnlel  is  useil.  "I'issetitially,  the  dominant  aspects  of  the  dynatnics  are  included  in  the  moilel, 
an<l  one  accounts  for  the  many  neglected  effects  by  inlrodticing  additional  nncertainty  into  the 
tiiodel”  [19:25].  The  pseudonoise  is  added  to  each  filter.  The  difficulty  of  usitig  this  approach  is 
how  mtich  pseudonoi.se  to  add.  If  too  much  p.seiidonoise  is  added,  as  in  the  ca.se  where  Qj  is  made 
artificially  large  to  account  for  the  inadec|uacy  of  the  propagation  model,  then  the  propagation 
equations  may  not  remain  usable  for  online  calculations.  4'he  addition  of  e.vcessirc  drir  ing  noi.se  to 
the  model  will  mask  the  difference  between  a  good  and  had  pro|)ag;ition  model.  The  addition  of  too 
little  dynamics  driving  noi.se  will  give  an  unrealistically  “perfect”  perspective  of  the  dynamics  model 
and  can.se  the  filter  to  use  the  model’s  characteristics  above  any  real-world  measurement  inputs. 
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'I'iii'  MMAE  algorilluu  lias  the  additional  hurden  of  using  this  same  dynamics  noise  tlironglioni 
the'  filter  s|)ace  for  the  impletnentation  used  in  this  research.  'I'he  tuning  of  Q,/  may  be  diliiciili  to 
achieve  with  a  single- valued  matrix. 

The  selection  of  adequate  measurement  noise  is  not  new  to  the  M .M.\ lv/.\l .M.\C  algcn  it  Inns, 
hut  is  a  Kalman  filter  problem  in  general.  Consider  Kqutition  (l.b’i)  which  provides  the  ICdman 
filter  gain,  [f  R  were  set  to  too  large  a  value,  K  (or  at  least  many  of  its  elements)  would  tend  to 
go  to  zero,  causing  tlie  state  propagation  equation  to  stop  using  the  residual  in  changing  tJie  \aliies 
of  the  states.  Likewise,  the  .state  covariance  would  al.so  remain  relatively  nnchangc'd  due  to  a  \eiy 
small  Kalman  filter  gain  value.  The  filter  must  “realize"  tliat  imjierfect  ions  in  the  dynamics  mode! 
e.\i.st  without  being  made  useless  in  the  process,  too. 

If  a  particular  filter  is  nearest  to  the  true  (larameti'r  location,  the  residnal  from  th;it  iiller 
would  be  expected  to  be  smallest  iti  magnitude  (relative'  to  tin'  lilier-coinpuled  covariance  Ai  of 
L'qnation.s  (2.12)  and  (2.13))  of  the  active  bank  of  filters.  Eijnation  !2.r2)  would  then  jirovide 
the  largest  conditional  density  value  for  this  filter  and  when  ajqilicd  to  Kqmition  (2.9),  would 
ultimately  yield  the  highest  probability  to  the  “best”  filter.  The  two  equations  mentioni'd  show 
that  the  filter’s  performance  is  dependent  on  the  difTerences  l)et\veen  the  residuals  of  the  filters  in 
the  bank.  Thus,  it  is  e.s.sential  to  ensure,  by  the  manner  of  discretizing  the  parameter  space  and  by 
t  he  means  of  tuning  each  elemental  filter,  that  the  residuals  of  “good”  versus  "bad”  fillers  have  very 
distinguishing  characteristics.  The  value  of  would  begin  to  dominate  however  if  the  residual 
values  were  of  the  same  magnitude.  Under  this  condition,  the  pk’s  are  dominated  by  the  since 
the  residual  values  are  equal.  The  conditional  density  then  becomes  dominated  by  the  filter  willi 
the  lowest  Ak  and  diverges.  “Unfortunately,  no  rigorous  general  proofs  are  available  concernitig 
till'  asymptotic  iiroperties  of  tlie  hypothesi.s  conditional  probabilities”  [19:  Uhl], 

Previous  theses  prevented  filter  lock  by  fixing  a  lower  bound  on  the  pk  [ti,  8,  1 1,  l.h,  21].  The 
hyiiothesis  conditional  probabilities  are  monitored  and  when  one  or  more  prolnibilities  fall  below 
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the  throsliold,  tliey  are  set  to  the  minimum  and  the  entire  bank  rescaled  to  maintain  tlie  unity  sum 
of  I  he  probabilities. 

2.  i  Moving-Bnnlc  Algorithm  Development 

The  MMAE  fdter  algorithm  presented  in  the  previous  section  presents  an  immense  iDinpu- 
taiional  Iturden  if  implemented  in  a  full  bank  form,  i.e,,  a  Piter  for  every  diM  icie  point  in  tin- 
parameter  sitace.  To  lessen  the  computational  loatling,  the  techniipie  of  a  nu)\  ing  bank  of  fillers 
has  been  investigated  in  the  previous  theses.  Maybeck  and  Ilentz  performctl  the  original  research 
that  has  direct  application  to  this  thesis  research  [22],  Their  work  showed  that  the  full  bank  of 
fillers  could  be  replaced  by  a  moving  bank  of  filters  that  could  follow  the  changes  in  the  parameters 
of  interest  [15:27],  The  techniques  examined  for  moving  this  smaller  bank  were  briclly  discussed  in 
Chapter  1  and  are  presented  here  in  detail. 

2.3.1  Moving  the  Bank.  The  moving-bank  MMAE  is  a  subset  of  the  full-bank  MM.\K. 
Conceptually,  the  bank  is  centered  around  the  best  estimate  of  the  parameter,  d'he  knowledge  of 
the  parameter  may  be  uncertain  a  priori  or  the  parameter  location  may  change.  I'ho  decision  logic 
u:icd  to  move  the  bank  is  critical  in  these  r.  ,!istic  cases  [8,  11,  15,  22]. 

“2.3.1. 1  Residual  .Monitoring.  The  effect  of  the  residuals  on  the  calculation  of  the  prob¬ 
ability  density  was  seen  in  the  exponential  term  of  Equation  (2.12),  The  exponential  term  contains 
a  (piadratic  factor  of  the  residual  vector  and  At,  a  function  of  Hi-,  Ri,  and  the  covariance  of  the 
state  estimate  of  the  A:-th  filter.  Let  this  term  be  defined  cis  the  likelihood  quotient  Z,i(C): 


Lk{ti)  =  ri(«,)Aj,  '(<i)ri(<i)  (2.16) 

Considered  in  the  scalar  case,  the  likelihood  quotient  is  the  ratio  of  the  residual  squared  divided 
by  the  filter-computed  variance  for  the  residual  [22:92].  Consider  the  effect  of  the  parameter 
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Iii  iiig  outside  till'  area  of  the  filter.  I'lie  residual  valiie.s  woiihl  all  become  larj;e  and  drive  ihi-  A 
likelihood  quotients  high.  A  preset  limit  could  he  determitiecl  through  p<‘rrormaucc  evaluations 
.such  that,  when  the  .smallest  Lh-  exceeds  some  value  T,  an  ajipropriate  filter  movement  action 
could  be  determined.  Additionally,  the  filter  closest  to  the  true  parameter  should  have'  the  smalh'st 
likelihood  quotient  and  thus  provide  an  indication  of  where  to  move  the  bank. 

The  use  of  residual  monitoring  is  limited  however  to  situations  where  thi'  system  is  not  subji ci 
to  "single  large  samples  of  mea,surenient  noise”  [22:92].  Ivxamination  of  10(|uatiou  (2.1b)  sIkjws  that 
the  values  of  Aif/,)  would  all  rise  appreciably  in  the  fac('  of  a  suddi’ii  high  value  of  v(/,)  which 
affects  the  residual  directly  through  the  realized  iiM'a.surement  value. 

2.3.  L~  Paranuttr  Position  Monitoring.  Equation  (2  M)  gave  an  expre.ssicni  for  the 
estimated  parameter  location  based  on  the  K  discretized  param<'t('r  locations  ami  hypoihi'sis  cciii- 
dilional  probabilities.  The  estimated  parameter  location  is  u.sed  in  this  technique  for  determining 
where  the  center  of  the  bank  should  be.  If  the  parameter  location  begins  to  move  to  the  edge  rtf  the 
ctirretit  bank  or  to  tnove  to  some  predetermined  distance  from  the  curretit  center,  then  the  decision 
t(j  tnove  the  bank  could  be  made. 

Reviewing  Erpiation  (2.H)  shows  tliat  tiie  strength  of  this  technique  is  that  it  j-  li('s  on  ;i 
history  of  measurements  rather  than  just  on  the  single  current  measurement  .  However,  whih'  this 
technique  is  not  prone  to  the  measurement  noise  problem  of  residual  monitoring,  it  does  rerpiire 
a  firiori  knowledge  of  the  parameter  location  or  at  least  sufficient  knowledge  to  |)lace  the  range  ol 
the  bank  around  the  true  parameter  location.  For  the  case  where  the  parameter  is  not  initially 
known,  the  bank  may  be  set  at  a  very  coarse  discretization  that  covers  the  entire  admissible  range 
of  i)aramcter  values,  and  then  changed  to  a  finer  discretization  as  information  about  the  parameter 
is  fleveloped  by  the  estimator. 
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2.3.1.  i  I’aniiiK  li.  r  Fosilitni  inid  ‘‘Velority"  t..‘<liinal(  Monitoiiuy.  Ihi^  lriliiiii|ii''  i." 
an  .'xtcu::.iou  of  tlic  previous  paramola'r  position  nionitoring  IcclmiciiK'.  If  llm  i>:ir;iiiictiT  is  vaj)  in^ 
slowly,  the  position  and  "velocity”  can  both  be  tnoiiitored.  'I'lie  parameter  velocity  is  the  chaiii^r 
in  parameter  location  between  sample  periods  or  over  an  extended  [leriod  of  time.  I'lnt  velocity 
vector  could  be  used  to  estimate  where  the  parameter  will  be  at  the  next  sam|di'  period  If  the 
imedicted  location  is  outside  the  current  bank  of  filters  or  beyond  a  certain  criti'rion  distaiua'  from 
the  current  center,  then  a  move  of  the  bank  is  indicated  in  the  direction  of  the  velocity  via  tor 

2.3.1. J1  IFohahilily  .Moniloruiy.  Moniloriii}!;  the  hy  imt  hesji  conditional  iiroliahilil  i's 
generated  by  Kipiation  (2.9)  provides  insight  into  which  filter  is  located  nearest  the  true'  parami'ti  r 
loc.ation,  Usitig  a  preset  threshold,  the  bank  can  be  moved  in  the  direction  of  the  filter  providing  the 
most  correct  parameter  a*,.,  i.e.,  the  one  with  the  highest  in-  value.  The  bank  seeks  to  center  itselfon 
the  filter  which  is  based  on  the  most  correct  assumed  parameter  value.  ,\s  with  parameter  posit  iein 
mtmitoring,  proltafiihty  monitorifig  uses  a  titne  history  of  measnrenient.s  and  i.s  less  .su.scepi  ible  to 
radical  changes  dtie  to  individual  large  realizations  of  rneastirement  noise'. 

2..3.£  Contraction  and  Expansion  of  the  Dank.  The  size  of  the  bank  need  not  be  fixed  and 
the  bank  filters  need  not  be  at  adjacent  finely  di.scretized  points  in  the  parameter  space.  If  the 
active  filters  are  not  associated  with  adjacent  discrete  parameter  values,  the  bank  behaves  as  a 
coarsely  discretized  moving  bank.  A  coarsely  discretized  bank  was  illustrated  in  Figure  1.7.  I'he 
estimates  provided  by  such  a  bank  may  not  provide  as  accurate  estimates  as  a  bank  where  the  active 
filters  occupy  adjacent  discrete  parameter  locations.  It  does  offer  a  higher  probalulity  that  the  true 
parameter  location  wall  be  within  the  bank  however.  This  is  a  good  attribute  for  initialization 
of  t  he  moving  bank  since  initial  convergence  has  a  higher  probability  of  occurring.  Maybeck  and 
Hentz  found  that  an  initial  coarse  di.scretization  of  the  baitk  improved  |)aramrter  aci|ui.sition  [22], 

The  bank  may  reipiire  expansion  for  the  case  of  a  jum|)  change  of  the  parameter  location 
Using  residual  monitoring,  the  jump  change  is  detected  and  the  bank  ex|>anded.  I'he  likelihood 
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ratios  of’ the  inipleineiited  filtf'rs  are  ex|)ected  to  be  large,  exceodiag  soiiu' t lir<'sluild ,  using  resiiliial 
iiioiiitoriiig  [11:2'.)].  Once  the  hank  Inis  been  ex|)anded  ;uid  a  ni-vv  ptirtiim'l er  location  isliiiiali' 
niade.  the  bank  iinay  then  be  cotitracted  around  the  new  location  and  the  algontlun  |)roceeds  with 
new  I'stinuitions  ;uid  control  inputs  as  before.  When  the  bank  is  coarsely  discretized,  monitoring 
the-  iiaratiieter  estinialioti  error  conditional  covariance.  P„,  provides  an  indication  wheti  to  contract 
tin'  bank.  When  the  norm  of  this  matrix  falls  below  some  chosen  threshold,  '  tin'  bank  can  be 
contrticted  about  the  paratneter  estimate”  [22:02]. 


Karnick  propo.sed  an  alternative  method  in  which  the  probability  of  a  side  ci  tin' 
monitored  such  that 


i'side('i)  = 


Eside/,(g(^)) 
Ed  .sides 


bank  tt'a." 


(2.17) 


is  the  probability  associated  with  each  side  [11:27-29].  The  t  hreshold  logic  for  cont  r<aci  ing  and 
expanding  the  bank  are  also  giveti.  If  the  side  probability  falls  Ix'low  some  t  hreshold.  I  hen  t  In'  side 
is  contracted,  ‘'('.'onver.sely,  if  the  probability  a.s.sociated  with  :i  side'  risi's  above  some  i  liri  shohl. 
the  remaining  three  sides  arc  •moved  in’  ”  [11:29].  The  third  possibility  involves  evaliniting  all  four 
sides,  atid  wiien  the  sum  of  the  sides’  probabilities  fall  below  some  threshold,  t  In'  bank  is  cont  rat  ted. 


2.3.3  Initialization  of  Nfiv  Elemental  Fillers.  The  issue  of  initialization  of  new  elenu'iital 
filters,  after  the  bank  has  been  expanded  or  contracted,  ha,s  been  studied  by  llentz  [8:29-20]  tunl 
Karnick  [l]:29-32].  Tlte  technicjues  developed  proved  to  be  comjnitationally  inlensi'.  Sinct-  there 
a()|)ears  to  be  no  lo.ss  of  performance  from  the  eipial  redistribution  of  the  discarded  filters'  prob¬ 
abilities  among  tlie  new  filters  [15:36]  ,  this  will  be  the  method  used  in  this  thesis  effort.  For  tlie 
case  where  all  nine  filters  are  new,  the  probability  will  be  divided  ecjnally,  i.e.,  pe  =  ^  for  all  k. 
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J..j  Slocha.^lic  (’dtittollcr  Design 

Both  tlu'  tnoving-haiik  and  fixed-haiik  MMAp]  algoritlinis  can  he  iisi'd  wilh  several  sluehasl  ic 
controller  di'signs  [20:9-20],  [22:92],  [8:32--10].  Tlie  “assniiual  certainty  (Xinivah'iice  design  "  nieihod- 
oKigy  presented  hy  Maybeck  [20;24 1-215]  will  be  used.  .As  a  form  of  fe<'dback  controlh'r  synllii  sis 
tiH-liidque.  this  method  s('parates  the  controller  into  the  cascade  of  an  estimator  and  a  determinisi  ic 
optimal  control  fimction.  This  method  is  also  known  “as  the  forci'd  or  heuristic  certainty  e(|niv- 
ahiice  design  teclinique"  [20:2'll].  While  this  technique  lead.s  naturally  to  .M  .M  .\  h-liasi-d  control 
(using  the  X  of  Fi”ure  1.2  |)remnll  iplied  hy  a  singh'  gain  — G'(a)  evalnalod  on  ihe  basis  .  lie 
estimated  paramt'ter),  the  implementation  used  in  this  tlu'.si.s  is  .M.M.AC  coiiirol.  which  lu'imanly 
differs  in  the  implementation  of  tlie  controller:  the  .struct  nre  shown  in  Figure  1,1. 

The  estimator  for  this  thesi.s  is  th<*  moving-bank  .MM.AF.  'I'he  .s.vstems  d(.weloped  for  each 
[t.aramcter  location  are  assumed  to  be  linear  and  time  invariant  with  stationary  noise.  'J'lie  I'stimatc' 
of  the  state  i.s  propagated  with  tlte  previously  developed  eqriations  (from  Chaj>ter  1,  Fqtiatiotis 
(1.10)  to  (1.15),  atul  in  this  cha|Her).  Likewise,  the  .system  is  updtitetl  using  tin.'  update  i'(pi,ii iotis 
developed  previously, 

tOach  controller  developed  is  a  linear,  quadratic  cost,  fnll-stat('  feedback  optimal  determiidstic 
controller  based  upon  a  specific  assumed  parameter  value  of  aj..  'I'he  cost  weighting  matrices  in  the 
fpiadratic  cost  are  constant,  and  the  error  state  space  formulation  is  time  invariant  [18:29],  Flu' 
output  of  the  controller  is  desired  to  be  the  optimal  control  function,  u",  such  that  the  c|u.adr,'it  ic 
cost  function 

J  -  E  ^  [x^(f.)X(/i)x(ti)  +  y^(<.)U(t,)u(<i)]  +  ^x^(t,v-(-i)X;x(t,v4.i)|  (2.18) 

is  minimized  [20:10],  The  variables  in  Equation  (2.18)  are: 
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J 

X 


=  cost  function  to  be  ininiiniz<Hl 
=  state  weighting  matrix 
Xy  =  final  state  weighting  matrix 
U  =  control  weighting  matrix 
.V  =  number  of  sample  ()erio(ls  from  la  to  i \ 
f.v+i  =  final  time 

t  v  =  last  time  a  control  is  applied  and 

held  constant  over  the  next  sample  period 


If  the  assumption  of  a  linear  system  driven  by  white  Gaussian  noise  is  used,  with  a  iiuadralic 
cost  function  given  by  Eciuation  (2.18).  then  the  optimal  consiani-gain  discreii'  linear  fecilback 
control  law,  a.ssiimiiig  full-state  access,  is  given  by: 


u-(f.)  =  -G-(f.)x(f.)  (2.19) 

The  full  state  access  is  replace  via  russumed  certainty  eqttiv'alence  by  the  state  e.slimale  provided  b.v 
the  .MMAL'  elemetital  filter  as.sociated  with  the  controller,  each  Ixhng  based  on  a  specific  parameier 
value  a*,  'fhe  controller  gain  is  given  by  [20:16]: 

G;(f,,)=  [U«i)  +  Bj(f.)K,(fi+,)Brf(/,)]'‘  [Bj(f,)K,(<.+i)$(f,  +  i,ft)]  {‘2.20) 

■.',’hich  is  also  part  of  the  solution  of  the  backward  Riccati  equation  solving  for  K,  (<,)  [20:l.a-!G]: 

K,:(f,)  =  X(fi  )  +  $^  (f-i+l  ,  <,)Kc(ti+l  )$(<i+l  ,  f, ) 

-[$^(f,  +  ,,f.)K,«.+i)Brf(t.)][U(f.)  +  Bj(f.,)K,-(f,  +  ,)B,,]-‘  (2.21) 

x[Bj(f.)Kc(<i+i)$(<i+i,<i)] 
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as  suivtal  liackwanl  from  the  termitial  comlition: 


Kc(^i  +  i)  —  Xy  I’J. 

1  hi'  selix'tion  of  the  weighting  tmitrires  found  in  K<ination  (2. IS)  is  ini|ioitaiit  in  deti  rminiiig 
the  jierfort  nance  of  I  lie  font  roller.  By  choosing  an  appropriate  X.  t  he  import  ance  of  maint  aiiii  iig  t  hi¬ 
st  nicture's  rigid  body  (losition  may  be  stressed  algoritliniically.  The  tise  of  a  diagonal  matri.x  with 
large  values  emphasizes  the  importance  of  maintaining  the  states  near  zero.  I  he  input  weighting 
tn  :t.rix  must  also  be  chosi'n  with  care  since  it  influences  the  amount  of  control  energy  the  controllei- 
is  allowed  to  use  in  order  to  accom|)lish  the  |>urpose  of  p'^oper  state  control,  the  degree  that  the 
itnlividual  states  are  allowed  to  vary  during  the  sample  |>eriod.  'liming  of  the  control  and  state 
Weighting  tnatrices  "is  usually  rcciuired,  in  analogy  to  timing  of  Qwf/,  )  and  in  the  Kaiinan 

filter  [20:10-11], 

'rite  relationship  given  in  Ecination  (2.19)  is  also  the  solution  for  the  deti.-rininisi  ic  IXJ  opiiin.al 
'ontrol  prohh.'in  with  no  driving  noise.  Since  we  are  assuming  that  the  s>sleiii  is  iiuise-ciiriupied 
and  only  noise  corrupted  measurements  are  acce.ssible  versus  |ierfect  knowledge  of  the  si  ate-,  the 
state  used  in  Equation  (2.19)  is  replaced  by  the  state  estimate  x(t*)  provided  by  the  ('li'ineiital 
filter  based  on  the  parameter  value  cik-  As  seen  in  Figure  l.-l,  once  the  state  estimation  has  been 
passed  to  the  controllers  and  the  hypothesis  conditional  probabilities  computed,  the  control  inputs 
are  multiplied  by  their  appropriate  probabilities  and  tlien  summed  to  provide  a  control  input  to 
the  system. 

2..')  Slim  III  ary 

'file  algorithms  for  the  moving-bank  MMAE  and  MMAC  have  been  develojied  in  ibis  rha|iter. 
The  moving-bank  MMAE  should  provide  a  substantial  computational  savings  compared  to  the  full- 
bank  implementation.  This  computational  cost  savings  is  carried  through  to  iIk-  ronirolli  r,  sinn' 


IK>(  nil  llu‘  ('iMiu-nfal  cuiitroi  iniMit.s  uiaHl  to  l)0  ralcdlatfd.  A  primarv’  asstim(>f lod  (hrc-ii^li  llii 
iKa -  [opulent  was  tliat  tin'  system  of  interest  rould  l>e  de'soribed  in  a  linear  manner.  W'itlioiH  tin 
a>si  1  mpt i< ‘11 ,  the  j'>rohlem  would  I'xH'ome  extremely  diflirult  to  address  with  an  (Jii-lme  rout  n i] f  r 
[  h--  next  ciiapter  lievrlops  tlie  linear  model  used  in  this  and  pr<-vi()us  tle-s.-s. 


III.  Iloliilini]  I  iro-liiiii  I  .Mndil 

I  I  n  1 1  <  "i  II I  lion 

111''  ''‘iiM'i.l  >  I  a  larj;i'  llcxihlr  s|)ai'i‘  slriirliir<'  is  l.tiii'j,  invrr-l  ioanil  ii'iii'j,  a  lU"  1  .i\  lia>' 
iii'"l''l.  I  ill'  slrurliirr  Is  all.iwnl  1(1  laitali-  aboiil  ,a  piviit  iioinl.  llic  hi'.li.  I  la  '■iij.mal  : 

III'"!'!.  !  I  iv '  si  I  aa  1 1  . 1  by  l.yiicli  ami  baiula  [li;K!  lb],  ib.i's  iiol  liava  ibis  |.i\' 'i  in'j,  I'  aiiii'  .  !  '  ' 
i|>iiamu'  (■liarm'li'i  1st  irs  an-  !  lioruiifibly  uikIitsIooiI.  'rhc  iisi-  ul'llu-  si  inci  ms  ,ii|.i,(i,  |  ,,,  i  ba;, 
a!  i'  iWs  I'  ir  t  la;  st  ml\  ul  lioi  li  ri',;iil  hotly  rolat  ion  and  bend  ins  luodc  d  uiaiim-.s  ;  I  'c  17,.  I  la  r.  a  a  1 1  n  ^ 
iii'  drl  lias  not  be  l  li  ell, angl'd  sinci'  it  wtts  di'vidopt'd  by  Karnick  [  1  1  ..S'.l-a.s] ,  and  iisi  ■!  be  l,,i~li|.  i 
[  1  ■  I  1 1  -  I';'.)]  and  \  an  Der  W I'rkcii  l2'^l:  1  l-bl)].  .-(nd  it  is  shown  in  1-  ignn'  d.  I .  I  bis  ebapt cr  pr.  s.  iii  s  i  bo 
iiiat'  rial  froiii  the  ap|ii'opriatt>  ebapters  orilic  l.asi  lltri'e  llii'ses. 

i.  Moili'U:  .'iii  tiiiil  ()nl<:  r  mill  Shill  ‘s’/n/ci 

.\  mat  In  111, 'll  ieal  luoib  i  is  |■l■(pli|■(■d  lo  di'tnniiini'  I  hr  frasibilil y  of  i  In-  M  .\!  ,\  l  i/.\l  M  ,\(  '  c  iiil  n  ,1 
' '  'hii  i(|iii'  lor  t  hr  large  s.paer  s,  rnel  iirr.  1  hr  second  order  diirrrrnl  iai  ri  pi  a  I  ion  "wii  ieh  .goer  ri  is  the 
lb  xui'.'il  e  il'ral  ions  of  a  stnteHirr'' [1 7:;i]  is  given  by  [ll;.■^h].  [b'.  IT],  [21.  1  I); 

.M'Hns  C\t)  +  l<v(n  =  F(/)  =  F|(0  +  F.,(/j  bt.l  ) 

win  I'r  the  bit  litiiid  Side  is  the  siiiii  of  the  products  of  liia.ss  and  accrlrrat  ion,  daiaping  and  ecloeily, 
and  stillness  and  position,  'file  rigfit  hand  side  of  I'biuafion  (li.i)  is  actually  coiipio.srd  of  a  drtrr- 
II  nil  1st  if  eoiiiponent ,  Fun.  I},  ;i.nd  a  noise  coiiipoiienl,  Fv(l).  As  nolt'd  in  the  above  r.xplanat  ion.  t  he 
i-apital  b  tiers  on  the  left  side  of  Ffination  (d.l)  represent  sprcilic  striietnrr  charai  l rrist ics  'b  liiird 

.'Is 

.U  —  constant  tisii  mass  matri.x 

f  ■  =  constant  ti.x ii  ilatiiping  mat ri.\ 

ill 


,8in' 

■  Non-structuro!  Mosses  - 

•  Acce  1  erofneters  and 
Thrusters  Co' located 

Numbers  -  Nodes  Letters  =  Structurol  Elements 


Figure  3.1.  Annotated  Rotating  Two-Bay  Trns.s  Model 
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I\  =  constant  Hxii  si ilfiicss  matrix 


1  h"  n-diiiH'iisional  r(/ )  vector  in  Krinat ion  (3.1)  n'jire.Si'iits  tin.'  posit  ion  or  pliy>ical  coonlinai  i-.- 
of  I  lie  St  nict  lire. 

If  tlie  external  disturbances  are  a.ssnmed  to  In'  white  tianssian  noises,  tlien  Kqnation  (3  1)  is 
cliangi'd  to; 

Mr{t)  +  Cr{t)  +  h'v(/)  =  —bii(t)  —  (;w  Ci  2) 

win 're  [l."):18] 


vi(/)  =  r-ditni.'tisiotial  vector  actuator  ini)Uts 

i  =  Hxr  iiiatfix  identifying  position  and  vt'locity 

rt'latioitsliip  Ix'tween  actuators  and  controlled  variables 
w  ,s-ditttensiotial  vector  representing  the  dytiatnics 

drivitig  ttoise,  where  s  i.s  the  nntiiber  of  tioise  itiimts 
<!  —  nxs  matrix  identifying  po.sifion  and  relalionsliips 

between  the  dytiatnics  driving  noise  and  controlled  lariables 


If  the  states  of  the  systetii  are  taken  tc  be  comoo.sed  of  the  posit  ioti  atid  velocity  vectors,  then 
Etpiation  (3.2)  may  be  written  in  state  space  form  as  [15:49]; 


X  =  Fx  +  D\\  +  Gvf 


(3.3) 


winre  the  states  are  stochastic  iirociesses  since  they  arc  driven  liy  tioise  as  well  as  a  delermiiiisl  ic 
input,  d'lie  stale  v('ctoi  is  given  by  [15:49]: 


X 


r 


r 


2nx  1 


(3  4) 


'I'lie  uncertain  partmielers  enter  llie  matliematical  model  via  the  matrices  in  Ivpiation  (3.3).  'I'lie 
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sysit'in  matrices  are  given  by: 


/•’  = 


-.v-'A' X,. 


J  2  »  X  2  f  j 


/i  = 


On 


G  = 


On  X» 
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!/n  X  5 


One  of  the  assumptions  taken  in  tiiis  thesis  is  tliat  the  measurement  noise  enters  the  sysiinn 
at  tlie  same  location  as  the  actuators  (thrusters).  Tliis  assumption  causes  tlie  (j  ami  b  matrici's, 
ami  therefore  G'  and  B,  to  be  equal  and  the  dimensions  rand  s  are  ec|ual.  too.  d'he  locations  of  the 
actuators  are  shown  in  Figure  3.i  and  are  discussed  in  detail  later  in  this  chapter. 

d'he  St  at<  s  are  posit  ion  and  velocity  contributions  of  the  various  rigid  hod\’  ami  bending  modes 
to  physical  position  and  velocity  variables  at  arbitrary  points  on  the  structure  which  are  measured 
by  accelerometers.  'I'he  accelerometer  outputs  are  integrated  once  for  velocity  and  integrated  again 
for  position.  .Since  the  velocity  and  position  measurements  are  derived  in  this  manner,  the  position 
and  velocity  measurements  are  co-located.  The  co-location  of  these  measurements  would  imidy  a 
highly  correlated  noise  measurement  matrix,  R,  and  corresponding  non-zero  off-diagonal  elements. 
However,  for  pur[)Oscs  of  simplicity,  R  will  be  assumed  to  be  diagonal.  The  actual  measurement  of 
tln  si'  v.dni's  would  be  transformed  via  an  analog-to-digital  conversion  (ADC)  so  that  the  coininiler- 
implemented  control  algorithm  can  be  used.  The  ADC  process  has  inherent  errors  due  to  finite 
wordlength.  For  the  [)nr|)oses  of  this  thesis,  the  ADC  error  is  assumed  to  be  .small  and  will  lie 
accounti'd  for  by  increasing  the  measurement  noise  matrix  R. 
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The  model  of  tlie  measurement  of  the  states  is  given  by  [ll:‘12],  [lh:50],  [‘i  h-IT]: 


where; 

(li  =  number  of  measurements 

y  =  an  uncertain  measurement  dislurl)ance  of  dimensioii 
III  and  modelcrl  as  a  zoro-m<>an  while  (iaiissian  nui.-Ne 
of  covariance  R 

Up  =  position  measurement  matrix 
//,  =  velocity  measurement  matrix 

Since  the  position  and  velocity  measurements  are  co-located,  Hp  —  //„.  'I'he  ecjuations  are 
highly  coupled,  that  is,  the  equations  are  not  independent.  This  situation  is  not  desirable  and 
therefore  modal  decomposition  is  used  to  transform  the  sets  of  c(|uations  into  iiHhqx.'ndent  modal 
equations. 

■3.  J  Modal  Decompositton 

The  equations  ju.st  developed  are  transformed  into  a  system  of  equations  that,  are  decoupled. 
‘■'I'he  general  response  of  a  complicated  system  can  be  broken  down  into  the;  sum  of  ii  simple' 
re'.si)on.se;.s’'  [3:2f50].  'The  assumption  that  the  matrices  are  constatil  allows  the  moeltil  de'romposit  imi 
to  be  useful.  If  the  matrices  were  variable,  then  the  advantage  of  using  modal  elexe)m|ieisit  ion  wouhl 
be  lost  [3:2(52]. 

Following  the  research  performed  by  Lynch  and  Banda,  the  damping  matrix  is  assumed  to 
be  a  linear  combination  of  the  ma.ss  and  stiffness  matrix  [Ifr-lj: 

C  =  aM^i3K  (3,9) 
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Due  to  the  transformation  from  physical  to  modal  coordinates,  the  actual  determination  of  o  and 
3  are  not  required. 

If  the  modal  coordinates  are  denoted  by  the  relationship  between  the  modal  and  physical 
coordinates  is  given  by  [17;5],  [11;42],  [15:51],  [24:48]: 


r  =  Tii 


where  T  is  tlie  n  x  n  transformation  matrix  composed  of  eigenvectors  determined  from  the  solution 
of  [1  l:43j. [15:51], [24:48],  [17:5]: 


i-MT  =  KT 


(■■hll) 


Note  that  K(|uatiou  (3.11)  relates  the  modal  frequency,  u-q  the  ma,ss  matrix,  and  the  stillness  matrix. 
The  damping  matrix  is  not  involved  in  the  solution  of  the  eigenvectors.  The  solution  of  a  and  J 
are  not  explicitly  found  since  they  are  not  required  for  the  modal  decomposition.  The  values  of  x 
that  satisfy  Equation  (3.11)  are  called  the  natural  or  modal  frequencies. 

Stibstituting  Equation  (3.10)  into  Equation  (3.3)  yields  the  transformed  slate  space  o(|uation 


to  [11:43],  [15:52],  [24:48]: 


X  =  Fx  +  Bu  +  Gw 


where  the  state  is  now  defined  eis  [11:43]: 


'I  lie  other  mairice.s  from  FCquation  (3.3)  are  also  transformed  and  calculated  as  [1 1:43]: 


(3.14) 


-T-‘M-*A'T  -T-‘Af-‘GT 

L  J  2rix2tl 
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B  = 


=  <7 


(3.15) 


The  open  loop  plant  matrix  may  also  be  written  in  terms  of  the  untlampccl  iiai  nr, il  rriMpiencic.s  ajid 
the  (lamping  ratio  of  the  t-th  mode.  /'  becomes  [25:1769]: 


F  = 


0  / 

—2Qu>i 


(3.16) 


Finally,  the  equation  describing  the  measurement  process  may  be  rewritten  as  fl  1:11]: 


7,  = 


HpT  0 

0  f/,T 


X  +  y 


(3.17) 


The  use  of  modal  coordinates  allows  the  following  assumptions  to  be  made  [17],  [11:11], 
[15:5‘1].  The  structural  damping  is  assumed  to  be  uniform  liirouglioui  the  siruciuia'.  I'he  dam|iiiig 
coefficients  selection  does  not  play  a  role  in  the  undamped  natural  frequencies  and  therefore  may 
be  selected  based  upon  design  requirements.  The  previous  tf.e.ses  used  a  valin'  of  C  =  (1.005  ha.sixi 
on  work  performed  by  Lynch  and  Banda  [17]  and  as  representative  of  many  space  structures. 

The  mathematics  of  the  model  have  been  developed  in  this  section.  The  parameters  used  in 
the  equations  of  this  section  are  derived  from  finite  element  analysis  of  the  structure.  The  next 
section  provides  a  detailed  explanation  of  the  structural  model  itself. 


3-4  Two- Bay  Truss 

'Fhe  two-bay  truss  has  been  discu.s.sed  in  general  terms  to  this  point.  'I'his  .si'clion  gi\c.s  a 
[diysiral  description  of  the  model  used  and  describes  the  sensors  and  their  locations.  Tin'  original 
truss  model  was  n.sed  to  study  the  effects  of  structural  optimization  on  optimal  control.  It  consisted 
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Table  3.1.  Structural  Member’s  Cross-Sectional  Areas 


Member 

Area  {in^) 

Member 

Area  {in~) 

a 

0.00321 

h 

0.00328 

b 

0.00100 

i 

0.00439 

c 

0.00321 

j 

0.00439 

d 

0.01049 

k 

0.20000 

e 

0.00100 

1 

0.20000 

f 

0.01049 

m 

0.20000 

g 

0.00328 

of  the  two-bay  truss  attached  to  a  fixed  body  [17:13]. 

■  l.-i-I  Background.  The  model  used  by  Lynch  and  Banda  was  inodifii.-d  liy  adding  noii- 
structural  masses  at  nodes  one  through  four  to  simulate  a  more  realistic  large  s|iace  structure'.  TIu' 
masses  model  various  types  of  loads  that  might  be  found  on  a  large  space  structure,  such  as  fuel 
tanks  and  control  jets.  The  model  was  further  modified  by  attaching  'he  structiiii'  to  a  hid)  witli 
a  mass  much  greater  than  the  truss.  This  configuration  approximates  a  large  space  platform  with 
appendages  extending  from  it.  This  configuration  also  allows  for  rigid  body  mot  ion  to  be  cxaminerl 

The  rigid  body  motion  is  established  by  letting  the  hub  node  remain  fixed  and  the  truss  move 
around  it  in  the  x-y  plane  as  shown  in  Figure  3.1  [11:47].  The  rods  connecting  the  truss  to  the  hub 
arc  much  larger  in  diameter  than  the  truss  members,  causing  this  link  to  be  very  stiff  compared  to 
the  truss  rigidity.  The  addition  of  this  physical  link  does  introduce  high  frecpiency  modes  into  the 
system  but  maintains  the  low  frequency  modes  of  the  original  fixed  structure  [15:.57]. 

3.4.2  Tivo-Bay  Truss  Construction.  The  thirteen  rods  that  make  up  the  truss  structure  arc 
listed  in  Table  3.1  with  their  cross-sectional  areas.  The  members  are  made  of  aluminum  which  has 
a  niodulus  of  ela.sticity  of  10  psi  and  weight  density  of  0.1  pounds  per  cubic  inch.  1  he  sizi's  of  ro<ls 
k.  1.  and  111.  were  arbitrarily  selected  to  obtain  a  stiff  link  between  the  truss  and  the  hub. 

d'he  non-St  ructural  masses  are  indicated  in  Figure  3.1  and  are  masses  of  1.291  Ih  —  sec-.  fliis 
tnass  is  very  large  cotnpared  to  the  member  mas.ses,  which  achieves  the  low  freciuencies  associated 
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with  largo  si>ace  structures  [17:14],  'I'lie  actual  mass  valuc.s  were  selected  Ijasi  d  on  an  optimizing 
technique  that  inaintains  the  lowest  undamped  natural  frequency  (associated  with  the  lowest  mode) 
of  0.5  Hz  [26]. 

Once  the  masses  were  determined,  finite  element  analysis  was  performeil  to  determine  the 
mass  and  stiffness  matrices  describing  the  model  [20].  The  physical  cliaracieri.siics  of  the  imss 
members,  length,  cro.ss-sectional  area,  modulus  of  elasticity,  and  weight  density,  wme  e\alnaii  il  lo 
determine  mass  and  stiffness  matrices  with  dimensions  equal  to  the  number  cd'  degrees  of  in  •■d.jm 
associated  with  the  model  [15:58],  Since  node  7  is  fixed,  the  elements  of  I  he  mat  rices  associated  n  it  h 
this  node  are  eliminated.  With  three  degrees  of  freedom  and  6  nodes,  the  full  inorlel  i.s  represc'iited 
by  IS  modes  and  36  states.  However,  only  motion  in  a  plane  i.s  being  considered  so  that  the  system 
may  be  fully  represented  by  12  modes  and  24  states. 

Seufiors  and  Actuators.  combination  of  gyros  and  accelerometers  is  used  to  provide 
motion  iiifornialion  to  the  control  system.  Accelerometers  are  located  at  nodes  1  and  2  as  shown 
in  Figure  3.1.  The  accelerometers  are  not  located  at  the  node  of  the  bending  morles  being  detected 
since  the  “displacement  of  the  truss  caused  by  the  bending  modes  can  not  he  deli'cied"  [15:50] 
under  those  conditions.  As  previously  stated,  the  output  of  the  accelerometers  are  integrated  onci' 
to  obtain  velocity  data  and  once  again  to  obtain  position  data.  Gyros  are  |d<ac('d  at,  tlu'  hull,  node 
7.  to  provide  angular  displacement  and  angular  velocity  information,  .\cluators  are  co-loc,Ui'(l 
with  the  accelerometers  at  nodes  1  and  2.  The  co-location  is  done  to  simplify  the  system  model, 
specifically  the  measurement  matrix  H.  The  actuators  at  nodes  1  and  2  are  assumed  to  be  thrusters. 
The  hub  al.so  contains  an  inertia  wheel  to  act  as  an  additional  actuator. 

3.4-4  Physical  System  Parameter  Uncertainty.  The  goal  of  this  thesis  is  to  determini'  whether 
.M.M.AC  can  control  a  physical  .system  with  uncertain  physical  parameters.  I  he  parameter  space 
is  composed  of  100  discrete  combinations  of  the  parameters  of  interest.  The  non-st  ruclural  ma.sses 
and  stiffness  matrix  are  varied.  'Phe  non-structiiral  mass  changes  are  motivated  by  relating  this 
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action  to  tlie  depletion  of  fuel  tanks  on  the  structure,  refueling  of  the  fuel  t.inks  (where  weight 
is  heitig  added),  or  weight  being  shifted  from  one  part  of  the  truss  to  another.  I'he  varitilioiis  in 
stiffness  of  the  truss  members  might  be  caused  by  structural  fatigue  in  a  rod  or  rods  or,  in  the  more 
extreme  case,  actual  failure  of  one  of  the  members. 

d.d  State  Reduction 

The  state  reduction  performed  on  the  system  of  equations  developed  follows  the  work  doni'  !>> 
Kokotovic,  ot.al.  The  need  to  reduce  the  state  dimension  is  due  to  the  large  ci.iiipiiia- 

tioual  load  that  a  21-state  system  would  place  on  the  MMAE/MMAC  algorithm.  .Additionally,  the 
actual  effect  of  the  higher  order  states  may  not  play  a  major  role  in  the  dynamics  of  the  system. 
This  issue  is  addressed  in  this  thesis  by  observing  the  effect  of  the  higher  order  slates  in  a  grttdual 
manner,  as  was  explained  in  Chapter  2.  The  truth  model  contains  all  2-1  states.  The  reduced  ordi'r 
model  is  developed  in  this  .section. 

3.5.1  Development.  The  general  deterministic  system  is  first  reformulated  as  [ll:52-.a7]: 

u  (TIS) 

where  xi  corresponds  to  “slower”  modes  to  be  maintained  in  the  design  model  ;md  xj  corresponds 
to  "faster”  modes  to  be  ignored,  and 


X] 

Fn 

F,2 

Xi 

X  = 

X2 

F21 

F22 

X2 

+ 

B2 

z 


Hx  H2 


X 


(3.19) 


If  steady  state  is  assumed  to  be  reached  instantaneously  by  the  “fast”  modes,  the  x-j  modes  are 
removed  while  maititaining  the  low  frequency  modes  in  Xj.  F\i  and  F22  are  square  tnatriccs  and 
Xo  can  be  expressed  in  terms  of  Xi  assuming  that  exists.  Then  the  higher  order  modes  may 
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Ill'  iwpre.sseJ  by  [2-1:58],  [15:(.;2],  [11:53]: 

=  F21X1  +  ^22X2  +  B2X1  =  0  (3.21)) 

X2  =  — + /?2«)  (3.21) 

Sulj.stitutirig  tliese  equation.?  into  Equations  (3.18)  and  (3.19)  yields  [23:58],  [15:62],  [11:53]: 

X)  =  (Fii  —  Ei2E22^^2I  )xi  +  (/?!  —  B2)u  (3  22) 

z  =  (//i  —  —  UnF-Tn^Bju  (3,2.3) 

The  term  multiplying  the  input  in  Eciuation  (3.23)  is  a  direct  feedforward  term  created  liy  the  m-der 
reduction  [15:63]. 

Applying  the  above  order  reduction  technique  to  the  original  truth  inodi'l  s}'>icii>  in  iiimhtl 
coordinates  yields  a  new  system  dynamics  matrix  given  by  [11:5-1],  [11:6-1],  [23:59]: 

0  /  ;  0  0 

— — 2^itvi  0  0 

.  (3.23) 

0  0:0/ 

0  0  :  — — 2(,'2^'2 

The  matrix  given  by  Equation  (3.24)  has  obvious  partitions.  The  upper  left  partition  contains  the 
low  frequeticy  modes  to  be  retained  while  the  lower  right  partition  contains  the  high  frequency 
modes  assumed  to  reach  steady  state  instantaneously.  The  partitions  correspond  to  the  F]i  and 
F'n  partitions  observed  in  Equation  (3.18).  Additionally,  the  submalrices  Fpj  and  /h)  are  ze  ro 
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Sulistituting  this  information  into  lujiiations  (3. IS)  and  (3.19)  yinlds  [15;61]: 

Xl  =  /’llXi  +  Bill  =  f'rXl  +  Bril  (3. '2')) 

7,  =  lllXi  —  H2F22  Ihll  =  HrXl  +  Orll  (3  2()) 

wlit're  the  subscript  r  denotes  “reduced-order.”  The  only  term  associateil  with  tin'  high  fretinency 
rnorles  is  the  direct  feedtlirough  term  in  Equation  (3.26).  The  other  matrices  an'  fomni  liy  trmn  aling 
the  states  associated  with  Xv  from  tlie  full  state  model. 

Of  the  terms  contained  in  Equations  (3.25)  and  (3.26),  Dr  is  the  most  coitiiilicatcd  since  tlie 
other  terms  are  readily  availalde.  Previoti.s  theses  have  shown  that  []l;55-56],  [15:65-66],  [2-};59-(.;0]: 


Dr 


-H2Fy/B2 


o 

1 _ 

[-w?][2C2U“.]  [-.e?]-' 

O 

0  H'r 

/  0 

// 

(3.27) 


whe're  “the  unmodeled  position  and  velocity  states  are  represented  by  He  and  //'  resjiect ively" 
[21:60]  The  terms  seen  in  Dr  show  that  it  is  dependent  on  the  state  terms  that  arc  a.ssnmed  to 
reacti  steady  state  immediately. 


3.5.2  Order  Reduction  Selection.  For  each  location  itt  the  discretized  iiaratiieter  space, 
eigenvalues  and  eigenvectors  of  the  unreduced  system  may  be  calculated  from  the  system  dytitimics 
matrix  F.  d'able  3.2  shows  the  average  nattiral  frequencies  and  damping  factors  associated  with 
each  of  the  12  modes  of  the  structure  for  the  100  filters  that  were  developed  for  the  discretiz<'d 
parameter  space.  'Hie  damping  factors  are  all  close  to  the  0.005  value  used  as  an  tipproximat ion  in 
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'['able  3.2.  Average  Damping  Factors  and  Natural  P'requencies 


.Mode  .Number 

Natural  Frequency  (llz) 

Damping  Factor 

1 

0.000000 

0.000000 

2 

1.413014 

0.007617 

3 

2.941387 

0.011020 

4 

4.953902 

0.005000 

5 

5.440887 

0.005790 

6 

7.078167 

0.004720 

7 

9.119647 

0.005000 

8 

9.872206 

0.005559 

9 

155.113838 

0.0051  IG 

10 

1418.481250 

0.005000 

11 

1811.221250 

0.005004 

12 

3138.732500 

0.00.5000 

pre\  iou.'i  thc.se.s.  'J  In.'  eigeiivaliie.s  fall  into  closely  spaced  groiijis.  For  ex:ini))l<',  oiu'  grou]i  of  iliree 
modes  is  seen  in  modes  4  througli  G.  The  thesis  by  Van  Der  VVerken  [21]  was  the  only  |ire\  i<5iis 
work  to  use  a  24-state  truth  model,  and  problems  arose  when  comparing  this  to  a  rediu.avl  ordi  r 
filter  motlel.  On  observation  of  the  frequencies  listed  in  Table  3.2,  moties  2  through  8  are  all  in  the 
same  decade  and  tln.'refore  may  all  play  a  role  in  the  performance  of  thi^  control  system,  altlioiigh 
the  amount  of  energy  in  each  mode  must  also  be  considered.  V’an  Der  VVerken  ele'ctr'd  tc;  usi'  the 
first  three  mod<.-.s  in  his  filter  model  corresponding  to  a  G-state  filti.'r  model  a,s  done  previously  Iw 
Karnick  [ll]  and  l,ashle(r  [1-5]  (versus  a  24-state  trulli  model).  As  explained  in  ('ha|>ter  1,  oin.'  of 
the  major  goals  of  this  thesis  is  to  investigate  this  mismatch  and  to  determine  the  nninl/er  of  states 
ne('(led  to  provide  ade(|uate  performance  from  the  MM.\F/MMA(.'  for  controlling  a  large  s|iaee 
structure. 

?.  tj  Sutnmnry 

The  21-state  truth  model  and  the  6-state  filter  model  that,  may  he  employed  in  a  MM  .VC  con¬ 
troller  have  been  developed  in  this  chapter.  Tlie  degree  that  the  filter  model  wiU  re(|nire  increased 
slates  to  [trovide  adequate  performance  is  the  purpo.se  of  this  thesis  and  tin'  simulation  plan  to 
determine  this  increase  is  presented  in  the  eext  chapter. 
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/!  .  Simuldlion 


I  luliddurlion 

1  lie  space  st  rile' t tin'  under  invest igat ion  is  simiilatf'd  in  order  to  st ndy  the  rehiist  nl'^s  ol'  l  In- 
I'st  iniatioii  tind  control  capiahilit  ii'S  of  the  MMAL'  aiui  MMAC  ti'rhni(]n('  eniplo_\'ed.  do  <h  li  [  [inie' 
the  statistical  [properties  of  tin'  est iiiiat ion/cont rol  proct'ss,  Monte  ('arlo  anaicsis  is  pe  il'oiii!'  I 
1  li'  ■  s,  )1 1  ware  iisi'd  has  lieen  di'velopi-d  ovi'r  a  [leriod  of  years  ami  is  <  ml  lims  1  i  ii  this  cha  [  .1  >  i  1  !.. 
siinulat  iciii  plan  is  presi'iited  at  the  einl  of  this  chapter. 

Mdiilf  C'liiiii  An<ilij!'is 

rin'  .Monte  Carlo  analysis  perfornied  hy  the  .siimilation  sofiwan'  pro\iihs  tlm  siatisijeal  in- 
lorination  ahoiit  the  perlbrinance  of  the  M.M.AH  and  MMAC  tdgi.irit  Inns  d>‘s(  rilipsl  in  the'  pn'\'i"iis 
chapti’i's.  If  the  system  tinder  invest igiit ion  wen'  fully  litii'ttr,  ihi'ii  ccAariaie'e  aiial'.sis  eeiild  Is 
nseil  [18:3'20].  Ilowover.  the  ada|)tive  nature  of  the  MMAE/MMAC  te<hni()Ues  used  to  control  the 
large  sjiace  structure  requires  the  employnient  of  Motite  Carlo  analysi.s  to  ohlain  nitiny  stiniples  of 
the  process  atid  evaluate  tin'  statistics  of  the  [trocess.  As  with  the  work  pi'rfortiied  (yv  \'an  Der 
\\VTken  [21],  the  sitiiiilat ion  involve.s  two  tiiodcl.s:  a  2d-.state  "truth  inoder'  for  acenratt'  siniiilal  i(  in 
and  [lerformnnci'  evaluation  and  ;i  fi-state  filter  model  for  algorithm  (h'sigii  [iiirposes.  1  h'-  lill'  i 
iiiiidel  is  the-  same  luie  used  hy  Ktirnick  [ll],  I-a.slih'e  [la],  and  Van  Her  \\’erkeii  [2  Ij  I  he  2  1-si:il'' 
truth  nil  idel  wa,.  (h've|o|ied  in  t  In'  last  thesis  hy  Van  l)er  Wi'rken  [21]. 

1'  igiire  (  1 . 1  )  shows  a  block  ditigram  of  the  Monte  (.'arlo  analysis  simulation  which  is  a  vai  lal  ii  m 
of  that  used  hy  \'an  f)er  W’erken  [2  l:ls'>-l>7].  d’he  variables  of  interest  are: 

•  x,(M:  the  truth  model  states;  2'1-dimensional  and  in  modal  coordinates,  represent  iii.g  tin 

rigid  body  mode  and  eh'veu  bending  modes 


■i-l 


•  x(/,):  i>t  iiii.Ut  s  of  till-  system  s!;Ues;  ()-ilirn<‘ii.si(jiial  and  in  nioilal  isn  adi  nates.  re|ii'eseiit  iiiy 
the  ri”id  luidy  and  the  fitst  ti\e  Ix-nding  inodes 

•  a, (/,.);  the  s.'ctor  represetit  itig  tin'  strncttire  mass  and  stifFtiess  paratneters 

•  !}(/,):  estitiiates  of  the  iinrertain  parameter  vector 

•  the  ('rror  iti  the  paratneti'r  estimate  (lefincd  as  =  a,(t,)  —  a{t,) 

•  tiie  error  in  the  system  estimate  defineil  in  Krtnation  (1-8)  as; 


4-2 


Figure  4.1.  System  Estimator  Simulation 
1  ^  L 

i:=  1 

Ihi;  siatislics  of  tliti  parameter  estimation  errors,  e^,  or  true  states,  x,  or  Hx,  (for  control  evnl- 
nations).  may  also  be  obtained  similarly  with  appropriate  substitutions  into  Equations  (d.2)  and 
(■1.3). 

While  previous  thesis  efforts  ([11],  [15],  [24])  viewed  the  estimation  of  tlic  states  as  the  error 
si;;iials.  the  goal  of  this  research  is  to  view  physically  meaningful  quantities.  Tlie  position  and 
vrli.rity  estimates  of  the  filter  bank  corresponding  to  nodes  1  and  2  and  the  hub  of  Figure'  3.1 
are  compare'd  to  tlie  truth-model-generated  position.s  and  velocities  (taken  to  be  “noiseless'’)  via 
till'  mampulatiou  of  tlie  truth  model  rnea-surernent  matrix.  The  manipulafion  i,s  of  the  form  of 
varying  A  to  oiiserve  tlie  effects  of  tlie  higher  order  bending  modes  of  tin'  striic!iir<'  on  the  alnlily 
of  til'-  estimation  ti  rhn.qin'  'o  provide  accuiate  position  and  velocity  eslimales.  Van  Di-r  Werken 
[2  1.i>V-i;8]  viewed  "he  truth  states  a.sso.  iated  with  the  rigid  borly  mode  and  t  he  first  twrr  binding 
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inoilfs  \i'rsus  tlic  lilter  states. 

I!ii\\e\er,  tlii.s  ill.  .^i.s  ;i.sks  a  iiiiTi'  a|>nropriate  <|ueslioti:  How  w.-ll  <loi  s  tin-  al^nritlini  i  siiiiiai.’ 
(aii'l  al.s.i  control)  lie'  true  .shape  of  the  sjrace  slnirture?  'I'lie  siiiiulatioii  softwar.'  was  mo.lilir.l  to 
alh.iw  the  structure  positions  and  velocities  at  nodes  1,  2,  and  d  to  lx-  calcnlaled  from  the  truth 
in..id(d  and  estimated  from  t  he  filter.  .A  vector  of  error  signals  is  then  determini'd  hy  siiht  ra.  l  in.y 
the  filter  estimates  from  tin  true  positions  and  velocities. 

1  he  ot  her  simulation  performed  is  the  same  a.s  l  he  above  o[)cn-loop  est  imat  or  simulal  icui .  I  ait 
wiili  ;i  controller  i midiuiienf ed  for  closed-loop  c'.st im.at ion  an<l  control,  d'he  coiilrol  simulainiu  is 
depi'.ti'd  in  figure'  (1.2).  The  error  signals  will  be  geni'rtti.ed  ;ind  evaluated  .statistically  as  for  tin' 
est  imtu  ion  st  inly.  .Again,  the  i.'rror  signals  are  the  act  nal  differences  of  the  posit  ions  of  the  st  t  mi  ure 
as  .|i  termined  from  the  truth  model  versus  the  estimated  position  provided  by  the  lilter  bank.  I  his 
's  si.gnifictmt  since  the  goal  of  the  controlh-r  is  to  quell  the  oscillations  of  tin'  structure  that  may  lie 
itninced  In'  mo\'ing  t  he  st  met  iiri',  changes  in  structural  charticlerist  ics,  or  chtinges  in  nnii-si  rnci  iiral 
masses  on  the  structure.  Therefore,  the  primary  (|uantilie.s  to  In'  st  at  isi  ir:dly  rNaluaied  ar.-  tin' 
ad  mil  positions  of  the  structure.  In  particular,  the  position  means  and  v.aritincei,  of  nodes  1  and  2 
will  be  det.i.’rmined  using  H,x, 

Sofhrarc 

4- 1. 1  Inlioduction.  T  he  software  ii.scd  in  this  research  was  sttirted  hy  Ileiitz  [S]  and  then 
modified  and  nseil  hy  Filios  [d],  Karnick  [11],  baslileo  [1.5].  and  V;in  l)er  Werken  [21].  The  work 
riorme.l  through  I, ashler  was  iierformed  mi  aCDCf'yber  computer  (a  non-.Al'  Fl  resonrec).  \'an 
l)i  r  \Nerken  [2'1:tj0-70]  moved  the  F’ORT'R.AN'  programs  to  an  fl.xsi  (fltrO  siiperminirompnier  (in 
house).  During  the  hegimiing  of  this  thesis  effort,  tlie  programs  were  moved  once  again  due  to 
re'.l  I  ii'-f  iiring  of  cominiler  ri'.sonrces.  'file  programs  nenv  reside  on  two  separate  conipnier  system^ 
within  .AI  l  I’:  galaxy,  an  KL.XSI  bdOO  siiperminicompnter,  and  blackbird,  a  DfC  \’.A \- 1 1 /7M.5 


Figure  4.2.  System  Controller  Simulation 


.sui*'“rminicornputer.  The  programs  used  are  a  preprocessor,  a  processor,  and  a  postprocessor. 

4-3.2  The  Preprocessor:  SETUPS. F.  The  preprocessor  generates  the  matrices  used  to  do- 
scriije  the  structure  over  a  range  of  ma,ss  and  stiffness  parameter  points.  The  matrices  generated  an: 
<5.  B,j,  Gti,  and  H  for  the  system  model;  and  Dr  given  by  Equation  (3.27);  the  Kalman  filter  gains; 
and  the  LQ  controller  gains  for  the  100-point  parameter  space.  An  input  file  is  used  to  input  the 
state  and  control  weighting  matrices,  the  meeisurement  noise  covariance  matri.x  (R),  the  dynamic 
driving  noise  matri.x  (Q),  and  the  time  variables  (start  time,  stop  time,  and  time  increment ).  The 
ma.ss  and  stiffness  ma*rices  that  describe  the  structure  are  also  contained  in  this  file.  The  parame¬ 
ter  space  is  determined  by  multiplying  both  the  mass  and  stiffness  matrices  by  ten  different  scalar 
values,  thus  generating  the  100  point  .space.  The  truth  models  and  the  redticed-state  filter  models 
ai'i'  contained  in  two  output  files  to  be  used  by  the  processor. 
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.j.  I.J  I'ht  I 'riir( s.sor:  GBSK.F.  'I'lio  priiiiary  prMCc.-sor  iisi's  tin'  'nl';  iii.ili'iii  lu  ;ali  .l  l.\ 
ill''  iJicpi'oci'ssor  111  pi-rfurni  tin'  niovinj^-liaiik  simulation  via  Moiitn  (_'arlo  analysis.  As  l  iiyni's  l  1 
aii'l  1.2  inilicatr,  the  pioei-ssor  proiia^ates  llio  true  sysli'in,  vvliioli  is  a  I'lill  2  1-staie  nnnlrl  nf  iln 
St  1  inture.  ami  uses  nieasiiroments  of  this  system  to  uiidate  tlie  estimator  (ami  coni  roller).  .At  the 
nieasiireineiit  input  lo  tin.'  estimator,  noise  is  introduced  in  the  tdrm  of  a  white  t  hiiissian  ii"isi 
\i  cior.  I'lie  sc.iftware  has  the  caiial.iilily  to  ]>crform  bank  expiansions,  coni ract ions,  and  movem.  ni.- 
according  to  tin.'  logic  described  in  Sc'clion  2.3.1.  'I'hese  functions  an'  used  once  l  he  stales  of  i  h' 
lilfer  hank  have  been  njidnti'd  and  various  internal  parameters  are  compian  d  to  pn  sei  i  hie.-hohls. 
.\l  the  end  (if  c'tich  I  iiiie  increment .  perl  ineni  dal  a  is  writ  ten  to  out  |in  t  tiles  to  he  post  pin- c  sm  d  . 

I  Ilf  I'lisIjiKh-cssuri  U  i'.S  V  Ll'.F .  Once  t  he  simnlat  ion  is  coinpleled ,  I  he  da  1  a  I  iiiisl  he 
reduced  and  put  inio  a  |ire.sental)le  format,  i.e.  grapdis  of  the  mean  ±  (nn,'  standard  dev  jaiioi;  of 
i|mini  it  ies  of  interest,  ^I'lie  most  important  values  oti.si.'rved  in  this  thesis  elforl  a.ie  ilm  ei-i'ois  in 
position  and  velocity  eslimalion  ;uid  the  c.orresiionding  true  positions  rnid  '.e|o,-iij,s  for  ■o-nir'il 
evahiatiim. 

Fill) iihition  I’hni 

Tin'  main  objective  of  this  thesis  efrort  is  lo  determine  the  robust ne.s.s  of  M.M.VIO  and  .AI.M.At  ’ 
If.)  the  ('(feels  of  the  modi'lling  niistnalch  caused  by  a  21-slnle  irulh  mcjdel  and  a  (hsiaie  liln  r 
niod'  l.  'Idle  dilliciillies  enconnler(.’d  liy  Van  I)er  Werken  after  he  iueri.a.ied  the  model  stales  lo  2  1 
ni  iy  hav('  lieen  rhie  to  the  mismalrh  In'twc'cii  the  2'1-sliile  trulh  incKlel  and  the  li-slaO'  Idler  mod'  I 
.A  means  louhserve  the  cap.'diility  of  I  he  MMAi;/.\IMAC  al  gorilhms  to  provid('  useful  informaiii.in 
fin  the  r.ise  (jf  the  esi  iniaUji)  and/or  control  (in  I  he  ca.se  of  tin'  coni  roller)  is  lo  inonitor  the  posil  ion 
ami  \'el'.icity  estimation  I'rrors  ami  iruc'  position  and  veh.icit i'-s  of  the  sirncinre.  d  hese  \aln'  s  will 
la  eni|jhasiz(  d,  1  he  simulations  performed  dealt  with  tin'  perforniance  of  the  estmialioii  ofih'' 
sinic.tnre  slia|ie  and  rigjil  liody  orientation,  and  the  perforiiiaiice  fift!!''  ronirol  fdgorilhpp  in  ijn' 


pr.sonce  of  a  liiglier  order  truth  model.  Additionally,  the  work  of  previous  flic.se.s  is  atteMij)l<'d  due 
to  I  lie  changes  ni.idc  to  the  software  atid  the  computer  systems  on  which  ilie  wiuk  is  |)erfunnrd.  in 
order  to  corrohortite  correct  iniplenn'titation  of  these  software  changes. 

The  Dither  Sigiiiil.  To  view  the  effects  of  tlie  difrer<'tit  tnodes  of  tin'  huge  s|iaci'  struc¬ 
ture  adetpiatoly,  a  dither  signal  is  applied  at  the  beginning  of  each  10  .secotid  Monte  Carlo  an.d\  sis 
run.  'fhe  dither  signal  o.xcites  the  system  model  and  enhances  parameter  identilictttioti  [l.');7  l].  Kar- 
tdek  [11:81]  and  Lashlee  stated  that  square-wave  dither  signal  was  used  with  tnagnit  iidis  ul 

lUO  ttml  dOO,  resitectively.  Vati  Der  WTrken  states  that  a  "dO  llz"  sipiare- wave  signal  of  inagnii  inh’ 
100  w;is  used  in  his  research.  In  reality,  the  software  was  written  to  provide  a  psendo-dO^  dither 
.signal  generated  by  it.sing  the  sign  of  sin(d0<)  times  the  dither  magtiitude,  which  creates  a  sipiata- 
w;ive  whose  polarity  is  determined  by  the  sine  function.  The  dither  created  iti  this  tnatitier.  with  the 
sample  period  of  0.05  secotids  in  u.se,  did  not  provide  a  symmetric  dither  sigtial  since  this  viohates 
.S'fiatmoti’s  sample  theorem.  The  huh  and  attached  structure  liatJ  at)  initial  (jiost-dii  her)  \el(icit\ 
instead  of  retnaitting  at  the  tnill  or  fire-dither  position. 

From  the  Nyrptist  criterion  ['2d:dd4-d3(5],  the  maximum  frerpn'ticy  of  the  dither  must  be 


where /is  the  highest  frequency  contaitied  in  the  signal  of  ititerest.  I’hc  saiitftle  [leriod  used  in  this 
and  luevioiis  thesis  research  is  0.05  seconds  which  corresponds  to  a  maximutn  sample  fre'qitency  of: 


1 

2  X  0.05 


=  10//C 


rad 
(5d - 


sec 


(1.5) 


From  Fable  (d.2),  the  fust  eight  avertige  frequencies  are  all  below  this  10  Hz  fta'ipiency  and  slnndd 
he  aderpjately  excited  Ity  the  dither  signal  u.scd  for  this  re.search.  The  srpiare  wave  u.sed  does  cont.iin 


hartiiotiics  of  the  futidamenta/  fre<|ue/icy  with  siilficient  energy  to  excite  the  higher  order  mode>  of 


tiu'  structure.  The  effect  of  the  higher  modes  of  (he  structure  may  well  be  negligible,  but  maj(jr 
software  modifications  would  be  re<niired  to  change  this  limitation. 


4-4-~  MM  A  E  Study.  Van  Dcr  Werken  attempted  to  observe  the  effects  of  higher  order  modes 
by  com|)aring  the  filter  performance  to  the  truth  model  via  a  transformation  matri.x.  d'liis  tin. sis 
effort  follows  his  recommendations  to  scale  the  partition  of  the  truth  model  tiietisurement  matrix. 
H,.  associated  with  the  higher  order  mode  states,  to  be  scaled  by  a  constant.  A.  The  value  of 
A  will  be  increased  gradually  from  zero  to  one,  ijnle.s.s  dra.stic  effects  are  observed.  This  gradual 
introduction  of  the  nnmodeh'd  statt's  into  the  ('stimalioii  process  will  allow  for  an  invest  ii;ai  nai  ..I 
whi  ther  tnore  states  need  to  b(>  added  to  the  estim.ator  to  providi-  rolnist  state  estimation  ol  the 
s|i.'ic(.'  structure. 

The  actual  fortnation  of  the  modified  H(  is  performed  by  forming  a  \  matrix  di'fitied  as 


A  = 


^6x6  OsxlS 
0lSx6  AIisxlS 


and  the  modified  H(  tnatrix  is  formed  by 
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H,  =H,A  (d.T) 

d'he  reason  .\  is  diagonally  partitioned  is  due  to  the  way  Van  Der  Wi'rki'ii  generated  tin'  .system 
matrices  in  the  pre|)rocessor  progratii  [2‘1:22d].  Uidike  the  expression  givi'ii  in  Ki|uation  (1.25).  tin- 
truth  states  were  modified  by  V'an  D<;r  Werken  so  that  the  first  six  states  repri'setit  the  rigid  body 
and  first  two  bending  modes  for  both  position  and  velocity.  'Jhiis,  the  implemented  H,  will  ap])e;ir 
as  iti  f’,<iuation  (l.dO),  but  with  the  two  inner  column  partitions  interchanged. 

Once  H,  has  been  formed,  the  measurement  of  the  true  states  will  be  modified  by  tin'  A 


value.  A  A  of  0  corresponds  to  a  truth  model  and  filter  mode  order  ttiatcli.  A  A  of  i  would  lie  the 


casr  wli.'ro  the  truth  model  is  a  full  2  l  state  model  and  thus  a  mismatch  hetweeii  model  and  (iller 
will  exist.  A  benchmark  is  run  with  the  truth  model  and  the  filter  hank  al  tin'  .same  /narameli'r 
location.  d'his  position  represent.s  the  nominal  position  in  the  filter  .sjiaci'  since  the  sc, alar 

multipliers  are  both  one  at  this  point  and  is  used  to  maintain  sonie  comparability  with  prerious 
thesis  elforts.  The  first  number  of  the  (rosilion  pair  is  the  mass  location  while  the  second  number 
represents  the  stiffness  location  in  the  filter  space  as  generated  by  the  preproci'ssor.  .Additionally, 
the  initial  probability  of  the  center  filter  is  set  to  1.  This  arrangement  gives  an  output  llial  is 
comparable  to  a  single  Kalman  filler  with  full  artificial  knowledge  about  tfie  irarameter  coiidiliuii 
.Additionally,  the  filter  is  not  allowed  to  move  by  this  arrangemeiii .  Once  a  heiichmark  has  hi  i  n 
eslahli.shed,  other  test.s  ar«'  performed  to  determine  the  effects  of  t  he  truth  model/filler  iiiismai  eh 

do  test  the  (.'stimation  (lerformance,  the  true  paramc'ler  i.s  fixeil  ai  repriseiil al ive  jio.siiidii 
(T.ii).  which  corresjioiids  to  the  mass  and  stiffness  matrices  being  iimiiodilied  while  geiieraling 
the  truth  model  system  matrices.  'I'he  center  of  th>^  filter  hank  is  placed  at  irosition  (a. -A).  This 
corresiioiids  to  a  mass  matrix  scaling  of  0.8  and  a  stiffness  matrix  scaling  of  0.0.  'flie  .\lonie  Carlo 
rims  are  performi'd.  with  varying  values  of  A,  with  the  filter  fixi'd  (a  worsi  nisi  analysis)  and  with 
the  filter  being  able  to  move'  within  the  biter  space.  The  ability  of  the  algorithm  to  move  the  lille'r 
hank  [iropcrly  is  essential  for  the  cases  of  varying  parameters  and  jump  changes  in  t  he  paranioters. 

Tf.d  MMAC  Study.  The  MMAC  study  is  performed  by  enabling  the  software  to  control 
the  structure  by  applying  control  forces.  As  with  the  MMAE  study,  the  goal  is  to  deteiiiiiiie  tin' 
degree  of  iiifluenre  of  the  higlier  stales  of  llie  truth  model  on  the  iierformaiice  of  the  O-slaP-  filtrr 
aiul  controller  model. 

The  controller  study  is  performed  with  both  fixed  and  moving  hank  coidigurat ions.  Addilioii- 
all>  .  the  cases  of  the  informed  filter  hank  (the  filter  bank  fixed  at  tlie  true  parameter  location)  and 
iniiiifoniiffl  filter  bank  (the  filter  bank  fixed  at  tlie  .some  parameti'r  location  oilier  tlian  the  inn' 
[larameter  location)  will  be  performed  to  determine  the  pi'rformaiice  cliaracleiist  ics  of  tlie  M.M.AC 


algorithm.  For  these  studies  and  the  MMAE  studies,  the  values  of  the  dynamics  noise  slrengiii 


Q  and  the  measurement  noise  covariance  R  will  remain  unchanged  for  the  various  elemental  lilter 
and  controller  models. 

4  -i  -i  Disiid  tidncc  Rcjeclioit.  Once  the  estimation  and  control  |)ro(i  ssi  s  have  hi'eii 
galed,  the  robustness  of  the  algorithms  is  be  determined  in  the  event  of  a  dist  urliance  oiici'  thr 
estimation  and  control  process  has  reached  an  intermediate  steady  state  condition  (as  determiin  d 
from  the  MMAC  study  results).  The  disturbance  is  be  of  the  form  of  a  dither  signal  applied  in  the 
same  manner  as  the  dither  signal  applied  at  the  beginning  of  eacli  simulation  run. 

4.0  Software  Modifications 

•Several  [>rol)lems  with  the  software  were  encountered  during  tlii.s  resc'arch  effort.  J'he  I.M.SI, 
FOli.'FR.W  subroutines  [10]  used  in  the  preprocessor  had  to  he  replaced  with  a  newer  version  of 
the  IMSL  Mathematics  Library.  An  N  +  1  problem  was  discovered  in  the  po.st processor.  Addi¬ 
tionally,  several  poorly  coded  .sections  of  various  subroutines  were  modified  to  be  readalih',  Thesi- 
modifications  are  now  described  in  .soiiu'  detail. 

Oiii..'  of  the  iirohlems  caii.sed  by  the  conversion  from  IMSL  Version  9.3  lo  l.MSL  Wrsion  10. (1 
is  (hat  the  calculation  of  the  Kalman  filter  gain  matrices  for  the  filter  spare  is  possibly  numerically 
unstable.  This  problem  asserts  itself  only  for  large  values  of  Q.  In  order  to  create  a  filler  s|)are 
with  confidence,  tiie  value  of  Q  was  reduced  by  an  order  of  magnitude  from  that  used  by  La.shli'i' 
[l.'):91].  'File  filter  matrices  were  generated  with 
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rill-  luunorical  instability  did  not  prevent  the  generation  of  iisefnl  results  but  does  deserve  some 
future  consideration. 

The  N+1  problem  discovered  in  the  postprocessor  did  not  adversely  affect  the  results  of 
tlil.s  thesis,  since  it  was  corrected,  and  most  likely  did  tiot  drastically  affect  the  previous  thesis 
I’fR.irts.  When  the  processor  writes  data  to  the  intermediate  files,  it  writes  at  every  time  incremi  iii. 
including  the  first  and  last  time  points:  0  and  10  seconds  respectively.  This  tirrangement  cri  .ues 
200  ititervals  since  the  time  incretnent  is  0.05  seconrls.  However,  the  number  of  tneasuremeni s  is 
2f)  1 .  The  .Monte  Carlo  analysis  t herefore  writes  201  metisn remen ts  for  each  run  of  the  analysis,  .mi  1 
tiot  200,  as  a.sstimod  by  the  previotis  thesis  effort  [2'Ti.  For  a  large  number  of  tmalysis  runs,  the  \4  1 
|irublem  just  described  could  be  significtuit,  especially  if  the  beginning  and  ending  magiiiiiides 
of  each  sitnulation  rtiti  are  greatly  different.  With  the  N  +  1  problem  still  in  the  software,  the 
post proce.ssor  read.s  in  200  measurements  for  each  analy.si.s.  'J’he  201st  dattt  point  is  then  read  iti 
as  the  first  data  poitit  in  the  next  analysis  run  and  so  on.  For  a  Monte  Carlo  analysis  performerl 
with  100  runs,  the  la,st  .set  of  measurements  read  from  the  intermediate  files  would  be  misindi'xed 
by  100  data  points  and  possibly  provide  totally  erroneous  insights  into  the  process  being  stiulii'd. 

Sunn})  ary 

Section  ‘1.2  explained  the  need  for  and  'he  use  of  Monte  Carlo  analysis  ;is  it  pertains  to 
this  thesis.  Sections  1,3.2  through  A.3A  rlescril^ed  how  the  simulation  of  MM.\F/MM.\C  of  a 
large  space  structure  is  performed.  Section  4A  explained  how  the  simulation  would  be  used  to 
investigate  the  mismatcli  between  the  truth  model  and  the  moving  bank  of  Kalman  lilters  or  LCJC 
cotit  rollers.  Finally,  Section  ‘1.5  presented  details  on  how  t  he  software  was  implemi’iitc'tl  .-uid  some  of 
the  [irobicms  enconiilcri'd  with  the  implementation.  The  next  chapter  rlescrilies  the  results  derivc-d 
from  this  simnlation. 
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'■.I  I II I  miluit  U'li 

1  he  results  di'tlie  Moiile  Carlo  analysis  arc  presenter!  in  this  chapter,  'I'he  intent  of  this  ihe>is 
is  to  show  tin;  rohusun'ss  of  the  inoving-hank  M  M  A  E/M  .MAC  algorithm  in  the  presence  of  iinmoil- 
eler!  states  ill  the  truth  niodi'i.  \'an  Der  Werken  ['i-lj  modified  the  software  lo  providr-  a  21  stale 
truth  model,  d'he  rlimension  of  the  filters  remained  at  O-staie.  This  re.setirch  deiiionstrtites  wln  ther 
or  not  that  six-state  design  model  can  yield  adequate  performance  as  the  (-(Feet  of  nnmoihded  higher 
ordi'r  true  bending  modes  is  gradually  changed  from  nonexistent  to  fully  presimt . 

■T.  2  Multiple  MoiUl  Aditpltre  F.stimalton  Sitidtj 

5.2.1  Diipliriil/oii  of  Pii.'il  Work'.  'Dte  dili)lication  of  p.ast  work  was  performed  in  a  soiiierr  hat 
siib jc’Ctivi.'  mariner.  The  results  Jichieva.'d  by  Van  Der  Wi'rken  were  not  repeatable  [21:S  l-2l)id,  I'lie 
softwari’  [iroblems  describerl  in  Si'ction  '1.5  may  have  causeil  tin,'  problems  his  woik  ilisphayed,  but 
there  is  no  diri'ct  r-viilence  of  this.  Lashlce  worked  with  a  (5-st:\t(.'  filter  and  truth  model  [iri;,Sll-Slj 
and  tlius  only  tlie  cases  of  A  =  0.0  are  useful  for  comparison,  ffased  on  t  fie  results  l.aslilee  presents, 
this  work  compares  very  favorably.  'I'lierefore,  the  intent  of  this  work  is  to  iicce|)t  this  diiplicat  ion 
(jf  l.ashlec’s  work  a.s  a  basis,  and  to  extend  beyond  it  to  determine  the  rolnistness  of  tlu'  .algorithms 
in  the  full  presence  of  higher  order  modes  contained  in  the  truth  model. 

In  the  jrrevious  theses,  the  difference  of  the  states  were  used  in  the  Monte  Carlo  ainil\sis. 
bashlee  used  an  error  signal  of  [15:72-7d]; 

<V(b)  =  -  x(/,)  (el) 

where  r  is  an  ii  x  ii,  matrix  "tC)  imike  I  he  dimensions  comp.alilile,"  for  I  Ills  t  he.sis,  t  he  erroi  signa Is 
will  b(  physically  sigiiificaiit  values  r.alher  than  the  modal  varia'des.  The  signals  of  ii.ierest  within 
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tlio  simulation  diagram,  I'igurr  l.l,  ;>re  the  noise-free  ineasiireinont:  of  tlie  positions  anil  veloiil  i's 
of  nodes  1,  2,  and  3.  The  |dots  that  follow  show  the  errors  in  position  and  velocity  est imat <-s  ul 
these  nodes,  as  well  as  [laranieter  estimates  and  filter  movement,  plot.s, 

.•),2.2  .l/.U.-l/-.  lie  ti  chin  ink',  do  judge  the  romi.stness  of  tlie  MMAIi  algorithm  adeipialily  , 

a  ln  nchm;irk  is  needed  to  compare  to  other  residts.  dhe  estimation  henchmark  is  gi  in  rati  ,!  I,y 
set  I  mg  A  I  o  zero  and  phicing  the  lilt  er  l>;\nk  and  |)arami-ler  ;it  locat  ion  (7,t>) ,  i.e  .  I  la-  mass  pa  r.amei ,  i 
is  located  at  Its  seM  iilli  discjeti'  Nalin-  out  of  ten  and  the  slilfness  m.atri.v  is  at  its  sixth  di'sr,  i, 
loc.ation  out  often,  d'he  henchmark  analysis  is  shown  in  f  igma-s  .ad  through  -ad.  .\ddit  ioii,all\ . 
the  prohahility  of  filter  five  in  the  bank  (the  center  of  the  finely  discretiz.i’d  bank)  is  set  to  one.  In 
(,'ssoncc,  this  is  the  case  of  an  artificially  informitd  (with  corri'ct  iiarameier  v.iliies),  non-;tda[Ui\(' 
filter,  ddie  (7,6)  location  corresponds  to  no  scaling  of  tlie  nia.ss  and  stiffne.ss  m.'ilrices,  i.e,  .'is 
deii'rmined  by  the  finite  element  aindysis.  The  10  Hz  rc|nare-w;!ve  dither  is  applied  for  the  first 
O.a  si'cotnls  in  order  to  excite  the  sysfeiii,  and  the  actnalor  output.s  are  ze-roed  alti'j'  this  period. 
Idle  dither  should  stimuhate  the  first  eight  modes  of  the  truth  model  (and  |)rob;ibly  stimnlaie  ilm 
ninth  through  twelfth  modes  adequately  as  well)  and  jirovide  an  aih-qiiati.'  means  of  performing 
parameter  and  state  estimation. 

fdgnrc  5.1  shows  the  ftenchmtirk  |)erforinance  of  the  arlifiriaily  informed,  iioa-adapi  i\e  lilier 
when  the  p.artuiieter  loc.atiiin  and  filter  loc.ation  match,  dhe  onh-r  ol  the  Iriilh  mod'l  and  lilier 
ni.aicli,  and  (inis  (lie  filter  should  (irovide  ojilimal  (lerform.ance.  .Alli'r  10  seconds,  the  erri  t  in 
estiiinite  of  till'  position  of  node  1  is  near  zero  and  tlie  staiuhird  deviation  is  ap|iro.\iin.ately  II.Dllii 
inches  or  O.OOG'/c,  of  the  structure  length.  The  error  in  the  velocity  estim.ate  of  node  1  is  coni|>.ar;i- 
tively  steady,  with  the  .signal  settling  after  the  dither  signal  i.s  removed.  Figure  5  2  shows  the  same 
mrormalion  for  node  2  with  similar  results.  Figure  5..d  shcjws  the  errors  in  the  position  and  vidocil y 
I'siimalion  for  node  .‘i.  the  structure  huh.  It  must  be  noted  that,  unlike  ihe  node  1  mnl  nede  2 
e.siimalion  errors,  the  nodi'  3  errors  are  in  radi.ans  and  radians  per  seennil.  Mine  the  hnli  is  fixed  in 


sp.u'i'  ami  only  roiatos.  Again,  the  small  error  signals  and  corn-spondingly  small  drviati'ins  w.iiild 
mdie.atr  that  tin'  estinnttion  process  is  provitling  good  infitrmat ion. 

5.J.1  (iii-AdapItie  Siiiglc-FilUr  Analysis  inib  Variable  X.  Figures  a.l  through  .a,!)  show 
tla-  ri'sults  of  the  A  variation  rims  under  the  same  comlitions  as  the  henchmark.  i.e.,  ftm  a  single 
artificially  informed  filter.  The  position  estimation  errors  for  nodes  1  and  2  are  essentialls  I  he 
same  for  A’s  of  0..a  and  1.0.  Som«'  high  frequency  effects  are  observable  hnt  .are  not  signific.-inl 
in  magnitude.  1  he  same  c.annot  he  s.aid  of  the  errors  in  the  velocity  estimates  of  noih’.'-  1  ami  2 
and  the  error  of  the  position  estimntt'  of  the  huh.  The  error  sigmals  dramaticall\  show  the  Ingh 
Ire.iui'iicy  effect  o|  admitting  some  of  the  higher  ordi'r  modes  into  the  measnri''’ienl  for  the  (iller, 
1  he  magnitude  of  t  he  errors  is  an  ordt'r  of  magnitude  higher  th.an  the  henchimirk,  which  ma.v  cause 
a  significant  |)rohlem  for  the  controller  if  the  tuning  process  is  not  properly  ]ierformed. 

S.i.Jf  .\oii-Ailni)ln'(  Single- Filter  If’orsf  Case  Analysis.  Figures  5.11)  through  5.12  provide 
a  wfvrst  cas('  analysis,  d  he  singh;  filter  is  set  to  tlie  (5.5)  pttrameter  pc.isition  and  is  made  to  he 
non-adaptive.  fhe  true'  ptirameter  is  located  at  (7.(5).  Tlie  same  dither  signal  is  used  to  e.vcitr 
the  system.  The  plots  show  a  dramatic  loss  in  estimation  ca|>ahility  with  the  deviations  being  an 
ordi'r  of  magnitude  larger  than  those  for  the  benchmark.  This  simulation  clearly  shows  (he  mt'd 
fm-  an  adaptive  filter:  mismatch  of  true  and  filter-a-ssiinied  parameter  valm's  yields  unacceptahle 
degradatirjti,  fhe  \alui'  of  A  was  akso  varied  in  this  worst  case  analysis.  'Fhe  ch.'inges  in  A  showed 
m>  ap[ireci;dile  changes  in  performance.  See  Figures  5.13  throiigh  5.15  for  the  ca-'C  ,,1  A  =  (),5 

1  he  pi.>.>ili.in  eslim  iti  error  data  from  the  above  simulations  ;ire  summari/'ed  in  fable  .5.1 
I  he  time  averages  ot  the  data  suggests  (hat  (he  higher  order  modes  observed  in  I  In'  non-/ei(,>  A 
c.'.'es  dci  not  significantly  affect  the  |ierformance  of  the  filter,  'fhe  worst,  c;ise  analvsis  did  show  llnit 
tin  re  IS  a  strong  need  for  adaptive  estimation  in  the  rase  of  mismatchc'd  filleieassiiiiied  and  tine 
parameters.  '1  In’  worst  case  analysis,  i  e..  the  ra.se  where  a  fi.xed  filler  is  not  located  at  the  line 

in  the  A  cfdiimn.  fhe  dal.a  piesenled  for  tin  wore' 
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parameter  location,  is  indicated  by  the  " 


case  analysis  was  g('nerateil  will)  A  =  0.0.  'I'lie  \.'locily  estiiiiatf'  errors,  shown  in  I'alile  '>/!.  did 
nol  provide  the  san]<'  iiisights  as  the  position  estiiinite  errors  since  tlie  veloi  iiy  .  stiniales  did  not 
change  consid(’rali|y  as  a  I'nnction  of  A,  I'lie  data  from  these  tables  are  also  presenied  in  lie-  b.ir 
ch.irts  of  Figures  -a.lii  through  -a. 21  in  order  to  enhance  coinpari.sons.  1  he  first  ihri  e  lignres  present 
the  data  for  the  cases  of  .M.M.AF  with  varying  A  while  the  la.st  three  pire’seni  the  same  data  \Mlh 
the  'worst  ca.se''  (mismatched  filter/parameter  an<i  A  =  0)  data  included  for  direct  com|)a’is(  in. 
rile  trends  in  the  pi'rformance  of  the  algorithms  in  llie  face  of  immodehMl  modes  is  ribvions  in  the 
velficity  estimate  errors.  .Additionally,  the  “'worst  case"  analysis  shows  a  definitr>  degrtidation  iti  the 
algorithm's  performance  due  to  mi.smatch  between  filter-a.ssnmed  aiul  real  parameter  cables, 

■'T.d.d  .\f<niiii/-li(nik-  MMAt.'  Aiialy.si.f.  The  ni'Xt  perfortitance  atialysis  of  tin'  M.M.AF  algo- 
ithtn  was  to  determine  the  ability  of  the  bank  to  move  withiti  the  filter  .s));ice  wheti  the  itiitial 
paratneter  estimate  atid  filter  location  differ  from  tlie  true  par.atiK'ter  location.  The  simitlat  ion  was 
inifi.'ifed  the  .satne  a.s  for  the  worst  ca.s'e  aiiaJy.si.s  except  t));it  tlie  bank  was  allow'i'd  to  move,  'fhe 
probalnlit  ies  of  all  the  lilteis  iti  the  liank  wer<'  itiititilized  to  7,  tm<l  the  lower  limit  of  the  prof.-i- 
biliiii.'s  wtts  set  to  O.d.').  The  valiK'  of  O.O-a  was  used  after  It  w'tis  di'ti.'rtnitierl  to  provide  tie'  best 
iticicenietit  of  tlt(.'  batik  ttnder  the  ettrrent  softvvare  configuralioti.  Residual  monitoritig  wtts  list'd  for 
the  movement  logic  and  the  move  threshold  sot  to  0.2h,  again  determined  in  an  ad  hoc  maiiiu'r 
b.'ised  on  jrerforniance.  The  results  of  the.se  simulations  are  shown  in  Figures  .5,22  through  .a, .'ill. 

fhe  results  are  considertibly  different  from  those  .seen  with  no  movement  tdlowed,  ,'is  in  thi' 
w'tirst  case  analysis.  During  the  10  second  run  time,  the  position  estiimites  nndeigo  considerable 
changes,  fhe  jumps  seen  in  the  posit  ion  I'sl  imations  tire  corri'lated  with  the  actntil  position  ol  the 
sti'inture  as  e.xpecled.  'fhe  true  position  of  node  i  seen  in  F'igure  .5.2(5.  lit''  mid-st  rnctnre  point, 
shows  the  effi'ct  fif  the  mi.xtnre  of  high  and  low  frei|uenry  bending  modes.  .At  about  2.2.  I dl,  I'.d, 
7 a'.,  and  '.t.ti  seconds  of  the  run  lim<',  the  position  “smooths"  ;is  the  bending  mode  ('om|ioneiils  go 
out  of  ph.'ise  and  lend  to  <'ancel  each  other,  'fhe  changes  m  the  position  estiin.ile  error  plot  lor 


1  si  t'iii  to  foiTt'hite  to  tlii's<‘  points.  'I  he  cliangos  in  tlu'  nodi'  2  p  tsition  esl iinat ion  (.'rroi's  ai'' 
not  nsibhj  cormlat I'd  to  tlio  node  2  position.  As  the  structure  hends  ;ind  iiiovi  s,  and  the  positions 
aiT  iui-a.-:ii;  'd  ironi  tlie  truMi  model  via  r  =  H(X,.  tlie  new  ineasnrement  inhirmal  ion  is  passed  to 
the  til.er.'.  and  acted  u]ion.  'The  ch.angi's  in  the  liank  location  are  evident  in  I'  lgure  5.25.  i  he  mass 
par.'inn'ter  estimate  and  the  actual  shifting  of  the  hank  |ierforni<'(.l  well,  lingering  in  the  vicinil> 
1. 1  the  true  mass  parameter  location  (mass  =  T).  'I'he  stiffness  parameter  estimate  [lerformed  less 
Wi  ll,  with  the  stiffric-ss  estimate  and  location  remaining  near  the  initial  stiffness  location  (itnlial 
stiffne.ss  =  5  \<’rsus  the  true  stiffness  =  l5),  after  drifting  in  the  wrong  direction  initially.  Ag.'iiti  tie 
cli.inges  in  the  parttmeler  estitn.ates  atid  the  movement  of  the  b.atik  can  he  correhaled  somewhat  to 
cleinges  iti  the  tictu.al  positioti  of  the  structure  as  with  the  correlat ioti  of  tie'  postlioti  i-st  itna!  ton 
ert'i  It'S  discttssed  above.  Figttres  5.27  throtigh  5.20  show  thi*  results  of  t  his  atialysis  iierforiited  wit  h  A 
s('t  to  0.5.  'File  chatige  iti  A  revealed  the  same  types  of  respon.se  chatiges  ns  seen  with  the  art  ilici.illv' 
itiforitiei.l .  nott-adapt  ive  filter.  I'he  ri'sults  of  the  A  =  0.5  analysis  makes  this  behavior  loss  cletir 
sitn'i.'  the  (terform.'uict'  is  revi-rsed,  i.e..  the  stilTness  performs  be' ter  thatt  tin.'  ttiass  esi  ittiat  iott  .atid 
Itto\ettletlt  process. 

In  bolh  cases,  A  =  O.tl  atid  A  =  0.5,  the  visible  iialti.rn  of  the  position  t\nd  velocity  estim.aiion 
errors  of  nodes  I  and  2  were  not,  eviilent  in  node  2,  the  hub.  'Fhis  is  logictd  since  ihc-  hnlt  h.a,^  much 
gre.ater  iiittss  and  would  not  be  moved  or  otherwise  ciian,ge  its  posdion  after  the  dither  sign.al  is 
rein.  A’ed. 
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(5,5).  Position  and  (hi  Velocity  Estimation  Errors.  A  =  0.5. 
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■j.  Z  Multiple  Modi  I  AdnpliVf:  ConliolUr  Study 


it..!.!  Duplieation  of  Paett  Wink.  As  witli  llio  MMAK  sliidy.  llic  (liipliral  ion  o|  |)a>l  work 
was  iKTlorincii  in  a  soini'wlial  snlpji'.-iive  nianiior  and  llu'  saiim  liii.ll .it ;„;is  w.  ir  oiix  ivi  d.  1  li. 
s.uiic  an;d_\sis  was  porioriiii'd  willi  llin  MMAC  algorillini  as  with  tlin  MMAIi  .al”oiit  Inn  .and  lli. 
followiiiLr  sniisia't  ions  pnrsiait  tlic  MMAC  study  results. 

•a.  y.d  .\[.\I.\C  lienchiiiiirk.  ’[in'  .MM.\C  benclitiiark  is  deteriiiiiu'd  tin'  s.aiin'  w.iy  .as  t'.  i  lin 
M.M.Ma  In'iiclunark  and  is  ('(inivalent  t.o  a  single  artificially  infornn'd  I,C)(1  ('.nil  roll.  r.  I  li.'  nit.  i' 
and  p.'U'.'uin'tnr  Kn  at  ion  arn  set  to  (T.ti).  d  in'  filter  is  .'in  art  if i. dally  iiilon in  d .  ii.  'ii  - ada |n  i \ .  .  i.  a  n. -  ni 
of  I  In'  linn  iiig-liank.  I’hc  results  of  tin'  In'iiclunark  an'  shown  in  l-dgnies  .'i  li  1  i  In.  .nnh  .a,:’.,')  I  n;iii. 
■a. if.')  shows  the  position  of  nodes  1  aiiil  2  ;ind  contains  the  most  critical  iidoriiiai  ion  concaninn  tin 
.M.M.\(.'  algoritliiii’s  ability  to  provide  control  of  the  structure.  .-Mso  note  that  the  hub  actuator 
in|iut,  shown  as  the  rlarker  line  iti  Figure  h.M.  is  zero  after  the  initial  dither  during  tin'  first  Q..a 
seconds. 

Figure  .a.dl  shows  that  the  position  estimate  error  is  maintained  below  iKKri  inches  for  th.' 
|ieriod  of  the  Monte  (..'arlo  analysis  run.  This  is  an  imiarovement  by  a  factor  of  three  over  iln 
.M.M.AF  benchmark,  d  he  position  estitnate  at  node  2  is  similarly  improved,  d  lie  velocity  esiimate> 
and  the  node  3  [losition  estimates  do  not  show  much  change  in  their  |)erform;uice  ove'r  the  M.M.M', 
algorithm.  Figure  •5.34  shows  the  control  inputs  as  determined  by  the  M.M.AC  algoritlmi. 

d  he  control  signal  is  activated  at  the  O.-h  second  point  of  the  run.  d  he  initial  iiiagnii  inles 
are  large  lint  soon  ipiell  down  as  the  controller  gains  control  of  tin'  structure.  .\s  the  pri)ce.-,s 
cont  nines,  the  algorithm  causes  the  control  in|)ut.s  to  grow  suddenly,  d  his  activity  may  be,  in  part, 
a  phenomenon  that  is  called  bursting.  By  bursting  the  control  inputs  periodically,  tlii'  cniitrolh  r 
is  able  to  enhance  parameter  identification  and  thereby  improve  estimates  of  the  states  and  tin' 
positions  and  velocities,  dhe  control  input  magnitudes  remain  at  le.ss  than  -a  lbs  for  most  of  the 
control  run  as  show'ii  in  Figure  .5.34. 
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Siinilf -( 'lint idIIi  v  Aiialyr^ts  inlh  WiruihU  A.  .'l  i'iCi  l  A.  l'j 

>lni\v  till'  results  of  tin-  cotilrol  perrormaiiee  study  with  variations  in  A.  I  h.-  iiul\'  aiiiiifi  ialil'- 
(lili'  renee  liir  tin-  variatiuus  is  lu  the  initial  control  ac(|uisition  jilias''.  Ih  Iwii  u  It  A  , -11111  I  >.  i  i.u.l 
1  lie  iiit rocluct ion  of  the  higher  frecjui-ncies  ohviously  cause  the  control  iiiiuils  to  In-  coiisii|i  r;il)|\ 
largt  r,  thus  causing  a  transiriit  in  the  position  and  velocity  estiuiate  errors,  .Surprisingly,  tin-  sl.ili 
(shape)  control  fcir  the  case  with  A  =  1  is  superior  to  the  ca.se  where  A  =  I). 

■5.:?../  S on- Adaptuf  Single  Controller  ll'or.s/  Cane  Aniilgsis.  As  with  the  .\IM.-\F  stiidw  :i 
Worst  case  .\IMAC  analysis  is  perforiiual  with  the  hank  in  a  misinronueil.  iion-adapt  i\i-  iiioih  .  i.e.. 
the  hatik  at  posit  ion  (-a.a ),  t  he  parameter  ;it  |iosition  (T.ti).  atid  the  hank  lixi-d ,  I'igni'i  s  .A.  It;  1  h  r-  >iigli 
A, Alt  show  the  performtuice  for  this  configuration.  'I'he  .MM.\K  worst  case  analysis  showed  that  the 
est iiiiat ioti  process  suffeta'd  cousidertihle  degradat.ion  due  to  this  couligiiral  ion.  I  he  nou-adapi  i\i 
control  worst  case  analysis  shows  that  this  configuration  is  a  catastrophic  oik'.  I'he  control  input 
rises  iti  an  expotieutial  inatitior  frotti  the  initial  control  turn-on  transient.  I'he  position  ;ind  vi.-locity 
estimation  errors  grow  in  the  same  fashion.  If  this  were  to  be  allowed  to  cont  inue  tie-  st  met  u  re  would 
surely  self-destruct.  Adaptivity  is  thus  seen  to  be  critical  for  cotitrol  of  this  system.  .-\n  increase  in 
A  does  provide  some  ittiprovcrnenl.  in  performance  but  still  leads  to  probable  catastrophic  belun  ior. 
'I’he  restilts  of  the  benchmark,  A  variations,  and  worst  case  analysis  are  stimmarized  in  I'tibles 
and  5.-1.  Figures  5,5-1  through  5.59  present  this  data  graphically  with  the  first  three  illustrating  the 
('fleet  of  A  variations.  The  hist  three  figures  present  the  same  data  with  the  "worst  case”  iucludt'd 
and  show  temporal  average  results  over  the  simulation  time.  'I’he  meaningful  iuf’ormatioti  cotitaiiied 
in  the.se  latter  grajihs  is  the  tinstable  behavior  of  tiie  system  under  the  "worst  case"  mismatched 
conditions. 

•5.  7. .5  Moring-Bank  MMAC  Analysis.  'I'he  next  performance  analysis  of  the  .MM.\('  .dgi)- 
rithm  was  to  determine  the  ability  of  the  bank  to  move  within  the  filter  spaci-  wlu'ii  the  imiial 
parameter  estimation  and  filter  location  differ  from  the  true  parami'ter  locatiou,  and  to  investigate 
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its  caiiaoity  to  iiiaiiil ;iiii  stahl('  control  oI'iIk'  structure.  'I'ln'  siinnl.at ion  w.as  initiali  il  tlir  saiin'  a.~ 
H.ir  the  worst  ca.so  analysis,  ('.xcept  th.at  the  hank  was  .allowed  to  m<)ve.  'riic  prohahilil  ics  ol  all 
the  tillers  in  the  hank  were  initialized  to  i  and  the  low('r  limit  of  the  |)rohal.ilil  ies  was  set  to  II. 
Ihoidii.al  nionitorini^  w.as  nsc'd  for  the  movement  logic  and  the  mo\c  threshold  s,  i  t.i  II ’J.).  I  he 
results  of  t  hes('  simiihitions  .are  shown  in  Figures  5.(50  through  5.<S1.  'lhe.se  lignres  re|ireseni  lour 
sets  of  figures  for  frnir  difh'rent  A  values:  0.0,  0.25,  0.5,  and  1.0.  Tin'  vtdnes  ol  A  weia,-  chosen  .at 
;ui  ad  hoc  h.asis.  'I  he  value  of  0.25,  not  used  in  the  MMAE  study,  was  chosen  tr^  .assure  that  no 
performance  indicator  was  missed.  .After  performing  nniiK'rous  Monte  Carlo  runs,  the  de  cision  l(,i 
usi'  plots  such  a.s  Figure  5.35  to  show  only  the  mean  position  of  tin’  siiiiclnri-  was  inaile  sine.-  tin' 
adilitional  data  of  ±  standard  deviation  was  not  informative  to  the  (h'gree  that  neci  ssitateil  iln  ir 
presentation,  i.e,,  they  were  very  small. 

While  the  "worst  case"  study  and  the  moving-hank  study  hoth  started  the  filler  ;il  a  paiiim- 
eter  location  difforent  from  the  true  parameter  location,  the  ability  of  the  h;mk  to  move  :\nd  hreate 
tlie  true  parameter  proved  to  be  essential.  The  moving-bank  MM.AC  study  showed  that  tin.-  filter 
htuik  could  perform  a.s  well  as  the  artificially  informed  non-adapt iv<' singh’-lillr.'r  controller.  Con¬ 
versely,  the  non-adaptive  mismatched  filter  performed  cat<a.strophic;illy,  wit  h  l  he  si  rnci  nre  In'i'orniiig 
unsttddc. 

The  performance  seen  in  these  plots  indicates  that  there  is  no  dramatic  degr;nl;ition  in  .ulaii- 
tive  control  a.s  A  is  increased  from  zero  to  one.  Thus,  there  would  he  little  benefit  to  be  g.aiiied  by 
increasing  the  number  of  states  of  the  filter  model  beyond  the  current  si,x. 

5..(  MMAC  Dishntiaiicc  Rejection  Study 

With  the  (lata  just  presented  .available,  a  disturbance  rejection  study  was  p('rformi'd.  Fin 
intent  of  this  study  is  to  determine  whether  tlie  controller  can  Tnaintain  control  of  the  structure  in 
the  face  of  hoth  the  model  mismatch  and  an  unanticipated  ext<’rnal  disturhaiiri’.  For  this  study. 
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iioti'  that  A  =  0.0.  1  li('  aliility  to  reject  a  <listiirl>aiic('  is  essential  for  a  c(jiilrol  sysO'iii  o|H  rai  in,;; 
in  an  environment  pertnrheO  hy  both  man  ami  nature.  'I'lie  MM.\C  stmly  was  reworkrO  with  a 
modification  i,o  the  control  inputs.  At  tlie  ‘l-second  point  of  each  run.  tlie  dillu  r  was  turned  on 
again  for  0..a  si'coiuls  with  a  magnitude  of  .aO  versus  10  for  the  initial  dither.  ( 'h.^ed  liuip  rmitrol 
was  maintained,  as  oppose'd  to  no  control  inputs  being  applieii  during  the  inili.il  ll.a  second  dither 
period,  ['igiires  .a. S  I  thrcjugli  .a. 89  sliow  the  results  of  this  study.  'I'he  results  of  the  distiiibeil 
system  follow  a  jiattt'rn  similar  to  that  seen  with  the  pri'vious  M.M.AC  results.  1  lu'  velocities  of 
neales  1  and  2,  and  the  position  of  the  hub,  node  3,  are  greatly  afh'cted  by  any  (list  iirbaiici  ,s.  I,, 
they  the  initial  dither  or  a  disturbance  during  the  control  procicss.  d'he  control  inputs  iti  l  i'.;iti' 
•a. 87  show  that  the  structure  is  brought  under  control  in  less  than  O.-a  si'cc.uids  after  iln’  dist  nrliaiici’ 
is  rt'iiioved.  In  Figure  .a. 87(1.)),  the  controller  has  added  its  own  control  in|)nts  to  the  dist  nili.ince 
signal  and  remains  high  for  a  short  time  after  the  disturbance  ends.  .Afli'r  the  dist  iirbatici’  was 
applic'd.  the  controller  continued  to  perform  well  and  maintained  control  of  the  structure  although 
the  disturbance  has  increased  the  amount  of  “vibration”  about  the  zero  null  position  for  both 
nodes  1  and  2.  The  position  and  velocity  estimation  errors  returned  back  to  their  pri'-dist.nrbance 
levels  after  about  1.5  .seconds  after  the  disturbance  was  removed,  ddie  [ua'vious  results  for  thi.' 
moving-bank  .MM.AC  with  A  =  0  should  be  compared  to  the.se  results,  d'lie  controller  reji'cted  the 
disturbance  well,  and  this  provides  insights  into  its  ability  to  liandle  system  disturbances. 

5. .5  Number  of  Monte  Carlo  Analysts  Runs 

Several  Monte  Carlo  analysis  were  performed  in  which  the  number  of  runs  was  sigiiifiraiit l\ 
iiicrea.sed.  The  results  did  not  justify  the  added  expense  of  increased  computation  times.  Figure 
5.90  is  an  example  of  a  100-run  Monte  Carlo  analysis.  This  is  a  repeat  of  t  he  dist  urbance  reject  ion 
study  performed  in  the  previous  section.  While  the  plots  shown  in  F'igure  5.8-1  and  5.90  are  not 
identical,  the  general  characteristics  are  similar  and  provide  the  same  fundamental  information. 
The  velocity  estimation  error  does  show  that  tlie  post-disturbance  transient  is  .somewhat  longer 
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than  for  tin.'  100- run  Monti'  Carlo  analysis.  More  significant ly,  the  sl.rncl nre  jeinains  unilei  cuni  1 1 >1 . 
Ohviotisly,  after  lOU  rnns.  the  addit ional  data  smooths  the  means  and  varitinees.  'I'hi.s  coinp.irison 
of  10  and  lOO-rtiti  .Monte  Carlo  analyses  provides  corroboration  of  tlie  validity  of  the  results  givi'ii 
in  this  chaiiter  wliich  were  performed  witli  10-rnn  Monte  Carlo  analysis. 

5.6  SiniiiiKiri/ 

The  results  preesonted  in  this  chapter  covered  the  performance  of  two  nlgorithnis:  the  .M  .M  If 
.and  .\IM.\C  algorithms.  I'lie  estimation  hy  the  moving-bank  MMAE  of  the  |Misilions  .and  vi  h m  ii  i.  s 
of  the  nodes  of  interest  not  drtislically  affected  by  the  presence  of  higher  ordn-  mod'  s  of  iln' 
truth  model.  When  faced  with  the  hick  of  true  parameter  information,  the  singh-  ni)n-inf')rm'  i| 
filti.-r  [terformed  poorly.  The  problem  of  a  non-adaptive  fi.xed  (ilter/coni  rolhi-  was  eva  n  iii'iii'  '■li'.arly 
demonstrtUed  in  tin'  ca.se  of  the  M.\I.\C  study,  in  wliich  the  parameter  mismatch  caused  insi .-diility 
of  the  closed  loop  system. 

The  performance  of  the  MM  AC  to  estimate  the  positions  and  velocities  of  the  nodes  w;is 
improved  over  the  MMAE  performance.  Tlie  structure  was  also  effectively  brought  back  to  the  null 
position.  Tlie  moving-bank  MNEAC  preformed  nearly  as  well  as  the  benchmark  and  maintained 
adequate  control  of  the  structure.  The  disturbance  rejection  showed  tlial  the  controller  could 
maintain  control  of  the  structure  in  tlie  face  of  a  strong  short  duration  disttirhtuicc  signal. 


Figure  5.31.  Benchmark  Estimation  Errors  with  Single-LQG-Controller,  Node  1  -  True  Parameter 
and  Filter  at  Mass  =  7  and  Stiffness  =  6.  (a)  Position  and  (b)  Velocity  Estimation 
Errors.  A  =  0. 
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F  igiiro  5..i.3.  F3enclun<'irk  Estimation  Errors  with  Single-LQG-Controller,  Node  3  -  True  Parameter 
and  Filter  at  Mass  =  7  and  StifTness  =  6.  (a)  Position  and  (h)  Velocity  Estimation 
Errors.  A  —  0. 
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rosition 


f  igure  5.30.  Control  Inputs  for  MMAC  with  A  =  0.5,  (a)  Mid-struriurf’  actuator,  (h)  Einl  of 
structure  and  hub  of  structure  actuators. 
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41.  Estimation  Errors  with  Single-LQG-Controller, 
at  Mass  =  7  and  Stiffness  =  6.  (a)  Position  c 


A  =  1.0. 
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Figiir<j  5.13.  Estimation  Errors  with  Single-LQG-Controiler,  Node  3  -  True  Parameter  and  Filter 
at  Mass  =  7  and  Stiffness  =  6.  (a)  Position  and  (b)  Velocity  Estimation  Errors. 
A  -1.0. 
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Riii-e  5.14.  Control  Inputs  for  MMAC  with  A  =  1.0.  (a)  .Mid-structure  actuator,  (b)  End  of 
structure  and  hub  of  structure  actuators. 
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Figure  5.  IG.  “Worst  Case”  Estimation  Errors  with  Single-LQG-Controller,  Node  1  -  I'rue  Param¬ 
eter  and  Filter  Mismatch,  (a)  Position  and  (b)  Velocity  Estimation  Errors.  A  =  0. 
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Figure  5.55.  Position  Estimation  Means  and  Standard  Deviations  for  Single-LQG-Controllcr  Case 
with  “Worst  Case”  Analysis,  Node  1:  (a)  Position  Estimates  with  Worst  Case  Anal¬ 
ysis  and  (b)  Velocity  Estimates  with  Worst  Case  Analysis 
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Figure  5.65,  Mean  Structure  Positions 
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Figure  0.69.  Control  Inputs  for  the  MMAC  Bank  Movement  Analysis,  A  =  0.25.  (a)  Mid-structure 
actuator,  (b)  End  of  structure  and  hub  of  structure  actuators. 
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Figliro  5.S9.  Mean  Structure  Positions  for  Disturbed  System,  (a)  Node  1  Mean  Po.sition  and 
Variances  and  (b)  Node  2  Mean  Position  and  Variances. 


5-96 


17.  ('oiiciu.->ions  and  Rtcommcndut ions 


h  /  I  nil  (III  ltd  toll 

l  ilt'  robustiii'ss  aiialysisof  the  .MM.VF!  and  MMAC  algorithms,  for  tlie  strurtiire  ami 

rt.r  the  state  mismatch  between  the  truth  model  and  filters,  indicates  that  the  elfects  of  uime  hIi  IciI 
slates  do  not  cause  serious  degradation  of  the  algorithms'  pt'rformance.  'The  I'esuli.'  iii-e^eni,  ,|  m 
('liaptt'r  5  clearly  indicate  that  the  [iroposed  reduced  order  models  and  resuliing  lilt'.-r/l  i  iii  li  melel 
mismatch  could  be  tolerated  by  the  algorithms  and  that  more  research  should  be  |>erfoniie(l  to 
providi'  further  understanding  of  the  abilitii'sof  these  algorithms  to  control  a  hirge  spttce  structure. 

fj.J  Conclusions 

The  results  of  the  .\1.M.-\F,  stttdy  showed  that,  while  the  eslitnation  prce.-i'ss  of  tlu'  alg'Uiilim 
w;is  si.'tisitivc  to  tin.'  ttiismatch  between  the  filter  and  truth  models,  the  sensiti\  ity  wa.s  tiot  siiliiciciii 
to  justify  ittcreasitig  the  filter  tttodel  tiiode  content.  The  cotnputat ional  biirdi  ii  ciiciied  b_\  guiiig 
beyotid  the  curretit  (3-state  tnodel  would  be  dramatic  since  it  vvotild  appear  that  the  addition  ul 
se\eral  states  would  be  required  for  a  significant  improvement  in  performance.  .Mdil  ioiitdlv.  the 
Worst  ca.se  analysis  showed  that  the  performance  of  the  mismatched-paraiiK'ter,  nou<'ulni)l i\ e  lilier 
w;is  very  poor.  Ihe  ability  of  the  bank,  in  the  fine  discretization  mode,  to  mo\i'  within  the 
[larameter  space  was  good  enough  to  provide  equivalent  performance  to  tlu'  art  ifici, ally  informed, 
non-adaptive  filtitr  used  for  the  benchmark. 

Idle  results  of  the  MM.AC  study  were  even  more  conclusive  in  determining  that  an  iucrea.'-e 
in  filter  state  size  is  not  recpiired  when  compared  to  a  ‘2'1-state  truth  modi'l.  I'he  perloi mama' 
of  the  benchmark  showed  a  very  tight  control  over  the  space  structure.  I'he  worst  case  an.d\sis 
showed  (hat  the  mismtilch  of  the  a.ssumed  parameter  in  a  single'  hQCi  cnulroller  could  h  ad  to  a 
catastrophic  failure  of  the  .structure  it.self.  The  ability  of  the  bank  to  move  in  the  paraiiieii  r  spa<  e 
w.is  sufficient  to  provide  adeipiate  control  of  the  large  space  struct, ure,  despiti'  iinmoih'led  i  lfecis. 
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Ilii'  disturhaiici'  rejocliou  study  of  the  MMAC  algoril  Inn  sliowed  lliat  tlu'  conlrolli'r  could  o,-|m,iii1 
to  a  largo  short  duration  disturbance  and  maintain  control  over  the  structure. 

1  lie  luoving-hank  M.MAK/MMACd  algorithms  used  in  this  thesis  provide  gix.Kl  coiitrid  um  i 
the  two-hay  truss  structure.  There  remain  several  topics  of  research  that  shoiihl  he  invest  ig.iied  to 
feline  the'  |)erform;tnce  of  the  algorithms.  The  next  section  outlines  recomiiu'iidat  imis  that  should 
lie  performed  in  future  re.search. 

li  1  I‘CC!>III)VC  IlllilllOIIS 

Th  e  research  performed  in  this  tliesis  was  hampered  l.iy  a  numlx'r  of  fact  (.us.  The  scop.'  ..fihe 
rose, arch  pt'rformetl  was  not  sufficient  tocover  all  areas  of  inleo'st,  .and  research  shouhl  h.'  continued 
to  .address  sonic  of  the  concerns  raised  in  this  and  previous  thesis  efforts  fdsing  the  same  sp.ace 
sti'ucture,  filter/controller  models,  and  truth  model,  the  following  recommondal  ions  ,aro  made; 

1 .  1  he  values  of  the  dynamics  noise  strength  Q  and  the  measurement  noise  covai  iaina'  R  mat  tdci  s 
ill  the  Kalman  filter  g.ains  should  be  varied  and  the  effect  of  these  shotdd  be  iiivi  stig.ai.'.l  in 
detail.  Special  attention  should  be  paid  in  the  generation  of  new  filti.n-  and  iniih  m.jdel 
matrices  because  of  the  numerical  instability  problem  suggested  in  Chapter  5.  The  .^ol'twar.' 
modification.s  to  allow  for  timing  Q  and  R,  should  be  made  as  siigge.stxd  by  Can  Ih  r  W.  rkeii 
[2-l:210].  The  Q  and  R  should  be  scalable  in  the  sam-.'  manniT  as  the  mass  and  siiHiiess 
matrices  by  iii[)ntting  base  matrices  and  an  array  of  scalar  multipliers.  The  issue  of  U')ii-/i’i.i 
off-diagonal  terms  in  the  measurement  noi.se  matrix  should  also  be  invi'st igated  [2  l:21()j. 

2.  Verify  that  the  state  weighting  matrices  arc  tuned  for  the  liest  possilile  performance  of  each 
elemental  LQG  cotit roller  within  the  MMAC  synthesis. 

;i.  Th  e  issue  of  changing  parameters  should  be  researched  thoroughly.  The  parameter  should  b.' 
allowed  to  vary  slowly  and  perform  jump  changes,  such  as  would  be  caused  by  fuel  depl.-lion 
or  .structural  failure  of  a  truss  component,  respectively. 
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1.  'I'lie  p;iraiiifti'r  space  discretization  should  also  he  reevaluated,  since  the  fa|jal)ilit\  of  ilie 
lilter  bank  to  move  to  the  area  of  the  trn<'  parameter  was  not  optimal.  In  particu  ar.  piopi  r 
([lerhaps  coarser)  discretization  of  the  stiffness  parameter  sjiace  should  receive  attenlu  ii  in 
future  research  efforts. 

5.  The  use  of  some  variation  of  residual  monitoring  might  be  considered  for  performing  moviiiieiii 
of  the  bank.  One  possible  means  may  be  residual  magnitude  monitoring,  i.e,  moniioriim 
the  quantity  as  opposed  to  as  in  Kqiialicui  (l.’-’O)  .nid  .1 

accom[)!i.shed  in  this  research. 

b,  Th('  software  was  corrected  where  appropriate.  Afli^r  beisig  u.s.  d  and  mo.lilicd  for  lis'  in  ra 
tions,  the  software  is  in  need  of  being  overhauled  and  thoroughly  doi.'imienied  (in-c.il.  i.  In 
addition,  several  [lortions  of  the  code,  most  notably  the  dither  control,  has  fieen  rew  iiiieii  in 
an  ad-hoc,  manner  to  provide  useful  information  in  this  thesis,  Tlucse  warr.tnt  furtlier  stmly. 
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